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Current upper bounds of the neutron electric dipole moment constrain the physically observable 
quantum chromodynamic (QCD) vacuum angle \8\ < 10 -11 . Since QCD explains vast experimental 
data from the 100 MeV scale to the TeV scale, it is better to explain this smallness of \9\ in the 
QCD framework, which is the strong CP problem. Now, there exist two plausible solutions to this 
problem, one of which leads to the existence of the very light axion. The axion decay constant 
window, 10 GeV < F a < 10 1 GeV for a 0(1) initial misalignment angle #i, has been obtained 
by astrophysical and cosmological data. For F a > 10 12 GeV with 6i < 0(1), axions may constitute 
a significant fraction of dark matter of the universe. The supersymmetrized axion solution of the 
strong CP problem introduces its superpartner the axino which might have affected the universe 
evolution significantly. Here, we review the very light axion (theory, supersymmetrization, and 
models) with the most recent particle, astrophysical and cosmological data, and present prospects 
for its discovery. 
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I. OVERVIEW 



Strong interaction phenomena has revealed that the 
discrete symmetries of charge conjugation C, parity 

V and time reversal T are separately good symmetries 
of nature. Therefore, quantum chro mody namics(QCD) 
based on the gauge group SU(3) C |24l . ll87l Han, Nambu 
(1965), Bardeen, Fritszch, Gell-Mann (1972)] must re- 
spect any combinations of these discrete symmetries C, 

V and T to be accepted as the theory of strong inter- 
actions. Among these discrete symmetries, the CV sym- 
metry is not necessarily respected in QCD due to the 
nonzero QCD vacuum angle 9, an issue known as the 
"strong CP problem" . Since QCD is so successful phe- 
nomenologically, any possible solution to "the strong 
CP problem" is expected to be realized in nature. Cur- 
rently the most attr active solution lead s to the existence 
of a very light axion [HI HM HU Hzl Kim (1979), Shif- 
man, Vainstein, Zakharov (1980), Dine, Fischler, Sred- 
nicki (1981), Zhitnitskii(1981)L Searches for QCD ax- 
ions generated from the Sun [id |203| . Andriamonje et. 
al. (CERN Axion Search Telescope (CAST) Collabo- 
ration, 2007), Inoue et. al. (Tokyo Axion Helioscope 
Collaboration , 200 8)1 and remnant axions from the early 
universe [59|, 13181 . Rosenberg (Axion Dark Matter Ex- 
periment ADMX) Collaboration, 2004), Carosi (AD MX 
Collaboration, 2007)] are presently ongoing. 

T he sto ry of axions started with the QCD U(l) prob- 
lem [358l Weinberg (1975)] which is now understood, 



having been solved by the 't Hooft dcterminental inter- 
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action [346j,l348|, 't Hooft (1976,1986)]. The determinen- 
tal interaction is shown pictorially as the left diagram of 
Fig. Q] and the solution is shown as the shaded right dia- 
gram. The strong interaction makes the quark bilinears 
condense with the vacuum expectation value (VEV) of 
order v ~ 260 MeV. The phase of this interaction 9 origi- 
nates from t he Q CD vacuum angle which is known to be 
physical [55l . l207l Callan, Dashcn, Gross (1976), Jackiw, 
Rcbbi (1976)] and contributes to the neutron electric 
dipole moment (NEDM) by an order 9 times the neutron 
size, which is absurdly large. Peccei and Quinn(PQ) ob- 
served that there exists a way to make 9 just a phase by 
introducing a symmetry, now called U(1)pq, and hence 
physical a mplitudes d o not depend on 9 as in the massless 
quark case [29 4 12951 Peccei, Quinn (1977)]. In the stan- 
dard model (SM) , this phase is a pscudoscalar Goldstonc 
boson calle d the 'axi on' in the multitude of Higgs fields 
as noted in [360|, IH Weinberg (1978), Wilczek (1978)]. 
If the PQ idea is completed with Fig. [TJ this axion is 
exactly massless (but observable), and 9 would behave 
'unphysically' by the freedom of choosing an appropriate 
axion VEV, which was the original PQ idea. However, 
there exists subleading terms, proportional to one power 
of m q , which close the quark lines with the current quark 
mass instead of a condensation. Then, an axion poten- 
tial develops, and the axion becomes a pscudo-Goldstonc 
boson. The axion solution of the strong CV problem is 
cosmological in that the axion VEV chooses 9 = at the 
minimum of this axion potential. The currently allowed 
axion is very light and long lived. 

The properties of the axion (denoted as a) are chiefly 
given by its decay constant F a which sets the scale of non- 
renormalizable axion interactions through a/F a . Initial 
axion searches placed F a far above the electroweak scale 
and additional stringent bounds on F a were obtained 
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FIG. 1: The determinental interaction of light quarks. Chiral 
symmetry breaking introduces the anomalous 7/ mass term 
from the quark condensations. 
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Kim (1987)]. Axion astrophysics started by (l29l . ll30l , D 



cus, Kolb, Teplitz, W agoner (1978, 1980)], using earlier 
ideas from [111, HM K. Sato, H. Sato (1975), K. Sato 
(1978)], now gives a stringent lower bound on the de- 
cay const ant, F a > 0.5 x 10 9 GeV, from the study of 
SN1987A lollHHl Turner (1990), Raffclt (1990)]. With 
this large decay constant, the axion flux from the Sun is 
a small fraction of the solar neutrino flux, but still may 
be detectable by the CAST experiment and by the Tokyo 
helioscope. 

It is known that very light axions with F a in the 10 12 
GeV region (axion mass in the micro-eV range) can com- 
pose so me part of cold dark matter(CDM) in the un- 
verse Q, [Hi HqI Preskill, Wise, Wilczck (1983), Ab- 
bott, Sikivie (1983), Dine, Fischlcr (1983)]. The exact 
amount of the axion CDM depends on the initial ax- 
ion misalignment angle 6\ at the time of axion creation 
when the universe temperature was around the axion de- 
cay constant, T ~ F a . This observation puts the very 
light axion on the list of leading CDM candidate parti- 
cles. If indeed these cosmic axions compose a significant 
fraction of CDM in the universe, they may be detectable 
by collecti ng a xion-converted photons in t he cavity type 
detectors \$M Sikivie (1983)] as tried by [HI [l85|, De- 
Panfilis et al. (RBF Collaboration, 1987), Hagmann et 
al. (Univ. of Florida, 1990)] and now continuing at the 
AD MX experiment. 

Cosmology with CDM was the le ading can didate of 
the early universe in the 1980s [45|, |246| . |364| . Blumen- 
thal, Faber, Primack, Rees (1984), Kolb, Turner (1990), 
Weinberg (2008)]. Since then this view has given way to 
the new cosmolog y with the discovery of dark energy (DE) 
in 1998 [IMElE Perlmutter et al. (1998), Ricss et. al. 
(1998)]. The current view of the dominant components 
of the universe is f2cDM — 0.23 and Q\ ~ .73 with 
only a few percent consisting of baryons [340l Spergel 
et al. (Wilkinson Microwave Anisotropy Probe(WMAP) 
Collaboration, 2007)]. The most plausible dark matter 
(DM) candidates at present are the lightest supcrsym- 
metric(SUSY) particlc(LSP), the axion, the axino, and 
the gravitino. Here, we will review mostly on the axion 
and its CDM-related possibility. 

The need for DM was suggested as early as the 1930s 



[339|,l38l|, Zwicky (1933), Smith (1936)]. Since then, evi- 
dences of nonluminous DM in the universe has been accu- 
mulating: examples of which include flat galactic rotation 
curves, Chandra satellite photos, and gravitational lens- 
ing effects. If the galactic bulge is the dominant mass in 
the galaxy, the velocity v of a star located at r from the 
center should rota te w ith v ~ r" 1 / 2 . But the observed 
flat rotation curve [27l[ See, for example, McGaugh et al. 
(2007), and references therein] violates this expectation 
and implies an extended mass in the halo as p(r) ~ 1/r 2 . 
Also, the Chandra observation of X-ray and gravitational 
lensing ima ges i mplies this matter profile around the bul- 
let cluster |l06t Clowe et al. (Chandra Col labo ration, 
2006)]. Circular gravitational lensing images |208l . Jee et 
al. (2007)] also support the existence of DM. The DM 
density around us the Solar system is usually taken as 
p DM ^ (0.3 - 0.45) GcV/cm 3 . 

Current CDM candidates belong to either incoherent 
particles or coherent oscillations of spin-0 fields. From 
this disctinction, the bosonic collective motion such as 
the axion can be considered as CDM. The popular inco- 
herent CDM particles are the weakly interacting massive 
particles (WIMP) or decay products of WIMPs. A more 
frequently used independent distinction is thermal relics 
and nonthermal relics, but there is no strict correspon- 
dence relation between the incoherent coherent particles 
and the thermal-nonthcrmal relics. WIMPs are massive 
particles with weak interaction cross sections, first dis- 
cussed in terms of a heavy neutrin o, co rresponding to 
the RHS crossing point of Fig. [^a) |258l Lee, Weinberg 
(1977)]. The LHS crossing point corresponds to 10 eV 
neutrino [HI, [263, Marx, Szalay (1972), Cowsik, Mc- 
Clelland (1972)]. WIMPs, such as the LSP, arc thermal 
relics when their number density is determined by the 
freezeout temperature and are called nonthermal relics if 
their number density is determined by the other mech- 
anism such as by the decay of heavier relics [93, Choi, 
Kim, Lee, Seto (2008)]. In Fig. [^b), we sketch the 
axion energy density in terms of the axion mass. The 
shape is a flipped one from that of Fig. [Ha), because in 
the axion case the low mass and high mass regions con- 
tributes f2 a from different physics, one from the vacuum 
misalignment and the other from the hot thermal relics. 

In addition to the heavy neutrino, SUSY with R-parity 
conservation allows the LSP to be just such a WIMP par- 
ticle. The LSP interaction is "weak" since the interaction 
mediators (SUSY particles) are supposed to be in the 100 
GeV range. For a WIMP to be a successful CDM candi- 
date, usually the interaction cro ss section at the time of 
decoupling needs to be around [2461 . l340l . Kolb, Turner 
(1990), Spergel et al. (WMAP Collaboration, 2007)] 
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with n„,h 2 ~ 0.113 ±0.009. 



(1) 



This is roughly the cross-section for the LSP from low 
energy SUSY, which is the reason why the DM commu- 
nity is so interested in the WIMP LSP. Some superweakly 
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FIG. 2\ The Lee- Weinberg type plots for (a) the neutrino 

Q, v h 2 [2461 ] and (b) the axion Q. a h 2 , where h is the present 
Huble constant in units of 100kms -1 Mpc -1 . The dash line 
of (a) is for Q. v h 2 = 0.113. In (b), it corresponds to the 
hadronic axion. The blue and red dashlines correspond to 
the CDM and hot DM limits, respectively. 



interac ting p articles such as gravitinos, axinos, and wim- 
pzillas [1051 . Chung, Kolb, Riotto (1999)] can be CDM 
candidates as well, but their cross sections do not fall in 
the range of Eq. |T]). The CDM candidate particles are 
shown in the a- m t versus mass plan e in Fig. [3] with mi- 
nor modification from that of |319L Roszkowski (2004)]. 
The incoherent fermions, such as the neutrino and left 
ends of the bars of the axino and gravitino correspond 
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FIG. 3: Some proposed particles in the interaction cross sec- 
tion versus the correspondi ng pa rticle mass rrii plane. The 
skeleton is taken from Ref. [319l ]. The dashed curves repre- 
sent schematic shapes of fij versus the corresponding particle 
mass rrii. The small red square-box corresponds to the hot 
DM hadronic axion. Two small outside squares (in the cyan 
and blue colors) in the axion region are marked to show the 
plausible (GUT and CDM) axions, respectively. The abun- 
dances of heavy axino, gravitino and wimpzilla depend on 
how inflation ends. 



to the left crossing points of Fig. [2]Ja). The rest, except 
the axion, corresponds more or less to the right cross- 
ing points of Fig. [UJa) with the reheating after inflation 
considered if necessary. Currently, there are experimen- 
tal efforts to discover the LSP as predicted by SUSY 
models. Direct cosmological searches are also ongoing 
[fl HI l35l l4J HI. Ill l20al259L Jungman, Kamionkowski, 
Griest (1996), Bertone, Hooper, Silk (2005), Bernabei et 
al. (DAMA collaboration, 2003, 2008), Lee et al. (KIMS 
Collaboration, 2007), Angle et al. (XENON Collabora- 
tion, 2008), Behnkc et al. (COUPP Collaboration, 2008), 
Ahmed et al. (CDMS Collaboration, 2008)]. At the 
LHC, the probable LSP mass ranges will be looked for 
by the neutralino decay producing the LSP. 

It is known that density perturbations began growing 
much earlier than the time of recombination in order to 
become large enough to form galaxies in the young uni- 
verse. For galaxy formation, therefore, DM is needed 
since proton density perturbations could not grow before 
the recombination time, but DM perturbations could. 
With DM, the equality point of radiation and matter 
energy densities can occur much earlier than the recom- 
bination time sin ce D M is not prohibited in collapsing 
by Silk damping [HI, Silk (1968)]. If the WIMP mass 
and interaction cross section fall in the region allowed by 
Eq. O, it can be CDM. If the LSP is the only CDM 



component, then the LSP mass would give one number 
for the DM density, which may not be the case. Thus, 
even if the LSP is contributing to the CDM density, we 
may need the axion to account for the correct amount 
of CDM around us. This is possible in the anthropic 
scenario of very light axions because it is equally prob- 
able for the initial axion misal ignment angle 9\ to take 
any value between and 7r [3451 . Tegmark, Aguirre, Rees, 
Wilczek (2006)]. 

Here we review the axion, which is probably the 
most interesting Nambu-Goldstone boson jl80l . I28lll282t 
Nambu (1960), Nambu, Jona-Lasinio (1961), Goldstonc 
(1961)], as well as related issues. In Sect. [TT] we dis- 
cuss the strong CV problem and its plausible solutions. 
In Sect. IIIII wc review the most attractive solution giv- 
ing the very light axion and present the axion theory in 
terms of possible axion couplings defined with ci, C2, and 
C3 which will be used throughout this review. In Sect. 
IIVI we present axion astrophysics and cosmology. Here 
we present a new number for the cosmic axion abun- 
dance in view of the recent more accurate data on light 
quark masses. In Sect. [V]wc summarize the axion detec- 
tion ideas and the ongoing axion detection experiments. 
In Sect. IVII wc summarize the proposed very light axion 
models, including supcrstring axions. Finally in Sect. 
IVIII we briefly discuss cosmology with the axino, the ax- 
ion's superpartner. 

If the axion was observed, it would mark one of the 
most profound elementary particle discoveries because it 
would confirm experimentally the instanton-based argu- 
ments of QCD. In addition, if it were shown to be consis- 
tent with a cosmologically significant amount of axions, 
the CDM idea of the bosonic collective-motion would also 
be confirmed experimentally. If SUSY is correct and the 
axion is the solution to the strong CV problem, it must 
affect the evolution of the universe as well. 



II. THE STRONG CP PROBLEM AND 
SOLUTIONS 



T here are good reviews on the strong CV problem 
[Z3,[223,[29l, Kim (1987), Cheng (1988), Peccei (1989)], 
and here we outline a few key points. QCD with SU(3) C 
gluons is a confining gauge theory with thre e ligh t quarks 
below 1 GeV and A QCD = 380 ± 60 MeV pM Grootc, 
Korncr, Schilcher, Nasrallah (1998)]. The classical gluon 
field equations have the instanton solution [30, Belavin, 
Ployakov, Schwartz, Tyupkin (1975)], 



±G>iy which carries the integer Pontryagin index 



G u 



if(r)g 1 {x)d fi g{x), /(r) 
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where the gauge coupling is absorbed in the gauge field 
and g(x) is a pure gauge form with G M „ ex 1/r 4 for a 
large r and p is the instanton size. The (anti-) instanton 
solution satisfies the (anti-) sclfduality condition G^ v = 
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where G a » u = \e^ p<7 G a pa . The classical solutions with 
q = — oo, • • • , — 1, 0, +1, • • • , +oo, introduces a new real 
number 9 which parametrizes the \9) vacuum, 



\e) = 



oo 
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(4) 



Since ns are integers, in view of Eq. ([3]), 9 is a periodic 
variable with period 2ir. It has been known that 9 is 
an observable parameter [55|, 12071 . Callan, Dashen, Gross 
(1976), Jackiw, Rebbi (1976)]. In the 9 vacuum, we must 
consider the V and T (or CV) violating interaction 
parametrized by 9 = 9q + 9 woa k, 1 



C = 9{GG} = 



64tt 2 



*-* \j.v^-* per 



(5) 



where the curly bracket carries 1/327T 2 , 9q is the angle 
given above the electroweak scale and 6* W eak is the value 
introduced by the electroweak CV violation. This ob- 
servable 9 has led to the so-called strong CV problem 
from the upper bound on the NEDM. For QCD to be- 
come a correct theory, this CV violation by QCD must 
be sufficiently suppressed. 



A. Neutron electric dipole moment 

The interaction ([5]) is the anomaly term 0, [37|, Adlcr 
(19 69), B ell, Jackiw (1969)] which is the basis for solv- 
ing J348L 't Hooft (1986)] the old U(l) problem of QCD 
[358l Weinberg (1975)]. The important size of instan- 
tons for ph ysics is near the scale where QCD becomes 
strong. In [HI, 't Hooft (1976)], 't Hooft has shown 
that the determinental interaction of light quarks carries 
the same global symmetry as that of Eq. |(SJ), and it is 
customary to use this light quark determinental interac- 
tion rather than treating the gluon interaction ([5]). The 
early estimates of NEDM proportional to 9 from the de- 
terminental interaction are 2. 7 x 10~ 16 fl ecm [23|, Baluni 
(1979)] and 3.6 x 10" 16 9 ecm [HI, Crewther, Di Vecchia, 
Veneziano, Wittcn (1979, 1980(E))]. Other estimates in 
different methods arc 11 x 10~ 1 6 6> e cm [6a . Cea, Nar- 
dulli (1984)], 1.2 x 1O" 16 ecm [n| Schnitzer (1984)], 
3 x 10~ 16 6>_ecm 12791 Musakhanov. Israilov (1984)], and 
5.5 x 10- 16 9 ecm[212|, Kanaya, Kobayashi (1981]]. Com- 
prehensive reviews on the NEDM exist [ll8l . 1 3 lk Dar 
(2000), Pospelov, Ritz (2005)]. Recently, the NEDM has 



1 With the canonical normalization of the gauge field, the RHS of 
Eq. (0 is multiplied by g%. 
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FIG. 4: The insertion of the CV violation effect by VEVs 
of 7T° and r] in (a). They can be transferred to one vertex 
shown as a bullet in (b). With this bullet, the CV violation 
is present by a mismatch between the CV conserving RHS 
vertex and CV violating LHS vertex. 
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for the proton electrin dipole moment, i.e. d n 
Now, we estimate the NEDM as 
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where the CV violating scalar coupli ng q^ NN (the bullet 
of Fig. |SJa)) is estimated in Ref. [l!5l . Crewther, Di 
Vecchia, Veneziano, Witten (1980(E))] as 



g-KNN 



2(ms - mT,)m u md 
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U{m u + m d )(2m s - m u - m d ) 

(7) 
where we take Z = m u /md ~ 0.48, m d w 4.9 MeV and 
m s /m d ~ 20.1. From Eq. ggj of Subsec. UTTBl wc 
estimate the CV violating scalar coupling as 
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FIG. 5: Diagrams contributing to the NEDM with the bul- 
let representing the CV violation effect. Diagram (a) is the 
physically observable contribution. 



been estimated in the hard wall AdS/ QCD model with 
one extra dimension, 1.08 x 1O~ 16 ecm [l98l Hong, Kim, 
Siwach, Ycc (2007)]. 

The diagrams contributing to the NEDM are re- 
stricted. The neutron magnetic dipole moment arises 
at one loop in the chiral perturbation theory. If wc 
treat this neutron magnetic dipole moment operator 
Manom«cr M!/ nF™ as a vertex, tree diagrams do not con- 
tribute to the NEDM, because the magnetic moment 
term has the same chiral transformation property as that 
of the mass term and hence by redefining an external neu- 
tron field one can remove the phases in the neutron mass 
and in the dipole moment operator together. 

Let the U(l) chiral transformation of quarks in the 
broken phase be encoded in the neutron mass term as 
m n nLe l{a '^' /f ^'- a '^ 0/f ' +5/2) n R + h.c. (e la ' § instead of 
e 3ta e b ecause the baryon octet is spin-^). The vacuum 
expectation values(VEVs) of ir° and rf are calculated in 
Subsec. IIIIBI The CV violation is present by a mismatch 
between the CV conserving right-hand sidc(RHS) vertex 
and CV violating left-hand side(LHS) vertex as shown in 
Fig. 01b). The mass term of Fig. 0|b) and the neutron 
magnetic dipole moment term of Fig. [5]Jb) have the same 
chiral transformation property and the phases appearing 
there can be simultaneously removed by redefining tir, 
for example. However, the phase appearing in Fig. EJa) 
cannot be removed by this phase redefinition and this 
contribution is physically observable. Since Fig. [5{a) 
is the physically observable NEDM, for proton a similar 
argument leads to the same magnitude and opposite sign 



Note that Eqs. J7J) and (JSJ) give a factor of -10 diffcr- 
enc e. Th e existing calculations vary within a factor of 10 
[65l Il98j . These old calculations depend on the various 
approximation methods used, but non of these estimated 
the VEV of 7T°. For example, for Eq. © Eq. (11) 
of Ref. |l!5l Crewther, Di Vecchia, Veneziano, Witten 
(1979)] uses the SU(3) symmetric baryon octet coupling 
due to the CV violating interaction. On the other hand, 
for Eq. ([HI) the ground state vacuum of the mesonic fields 
has been used. After integrating out baryons, we look for 
the vacuum below the chiral symmetry scale. Then, the 
correct vacuum choice adds the value (J8J) to the value (7). 
But here we choose the one-order-larger value from the 
mesonic vacuum shift value (|5J) for an order of magnitude 
estimate, not to worry about the signs of the contribu- 
tions. So, we estimate the NEDM as 4.5 x 10~ 15 9 ecm 
fromEq. flSJ). 

Since the recent upper bound on the NEDM is \d n \ < 
2.9 x 10~ 26 ecm [H, Baker et. al. (2006)], we must require 



|0| <0.7x 10" 



(9) 



This extremely small upper bound on 9 has led to the 
so-called 'strong CV problem'. \0\ < 10 -11 is perfectly 
allowed but its tiny value is not explained given that it 
could have chosen a value anywhere between and — n. 
The strong CV problem is the quest to understand more 
satisfactorily, "Why 9 is so unnaturally small?" . 



B. Possible solutions 



In the remainder of this paper, we simplify the nota- 
tion, replacing 9 by 9 since there will not be much con- 
fusion. There are three explanations for the smallness of 



8 in the naturalness framework: 

Case 1. Calculable 8, 
Case 2. Massless up quark, 
Case 3. Axion. 

In this subsection, we discuss Cases 1 and 2 briefly, and 
concentrate on Case 3 in the subsequent sections. 



1. Calculable 6 

Naturalness of a the ory w ith a parameter (3 is elegantly 
defined by 't Hooft in [347|, 't Hooft (1979)]: The theory 
is natural if the symmetry of the theory increases in the 
limit of the vanishing /3. A frequently quoted example 
is the Dirac fermion mass, mipLipR + h.c, where m — » 
introduces a chiral symmetry -0 — > e l/375, in the theory. 

Regarding the strong CP problem, the appropriate 
symmetry is parity P or CP since the interaction ([5]) 
violates parity P, time reversal T and CP, but con- 
serves charge conjugation C. Requiring CP invariancc 
in the Lagrangian is equivalent to setting 8q at zero. 
However, the observed weak interaction phenomena ex- 
hibit the weak CP symmetry violations in the neutral 
K meson system and B — > K + n~ decay [9], Amsler et al. 
(2008)], and hence the needed introduction of CP viola- 
tion in weak interactions with 8q = must be achieved 
spontaneously. In this process one necessarily introduces 
#weak part in which can be calculated and required to 
be sufficiently small within the bou nd g ive n in (Q . Along 
this line, many ideas were proposed [27l l33 . l273 . l328l Beg, 
Tsao (1978), Mohapatra, Senjanovic (1978), Segre, Wel- 
don (1979), Barr, Langackcr (1979)]. This naturalness 
idea may be ex tended so as to only effect renormalizable 
couplings |l65l Georgi (1978)]. In any case, the introduc- 
tion of weak CP violation by spontaneous mec hanisms 
[26ll . Lee (1973)] or by soft scalar masses }l65l Georgi 
(1978)] must be checked with various weak phenomena. 
The current weak CP violation data fit s nic ely with the 
Kobayashi-Maskawa type CP violation J245I . Kobayashi, 
Maskawa (1973)], and these drastically different sponta- 
neous weak CP violation ideas are probably difficu lt to 
fit the data but are not considered ruled out yet [l93l 
He (2007)1 even though the spontaneous CP viol ation 
scheme [48|, Branco (1980)] in the Weinberg model [359l 
Weinberg (1976)] is ruled out [H, Chang, He, McKellar 
(2001)]. It should be noted, though, that the above pro- 
posed models have the difficulty in satisfying the bounds 



However, the Nelson-Barr type weak CP violation 
is mimics the Kobayashi-Maskawa type CP violation 
even though the fund amental reason for CP violation is 
spontaneous [H, [HI Nelson (1984), Barr (1984)]. The 
scheme is designed such that the Yukawa couplings are 
real, i.e. 0q = from the CV invariance. Next, the 
spontaneous CP violation is introduced through the sin- 
glet VEVs, which is the key difference from the previous 



calculable models. Thus, the spontaneous CP violation 
is required to occur much above the weak scale through 
the singlet VEVs, mediating it to light quarks through 
mixing with vector-like heavy quarks. In modern terms, 
the heavy quarks can be considered as the mediation sec- 
tor. Then, integrating out heavy fields, we obtain the SM 
quarks with the Kobayashi-Maskawa type weak CP vi- 
olation. To ensure Arg.Det.Mq — at tree level, the 
specific forms for the Higgs couplings to the SM quarks 
and the superheavy vectorlike quarks are needed. Be- 
yond the tree level, howeve r 8 is generated at one loop, 
typically with the form [50, 1 1781 . Bento, Branco, Parada 
(1991), Coffin, Segre, Welson (1980)], 



1 



''weak 
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-A/ 2 N^ (loop integrals) 



(10) 



where A/ 2 is the product of couplings and the Feynman 
loop integral is of 0(1). To satisfy the bound ©, the 
small coupling A/ 2 is needed. Some mechanism such as 
the famil y symm etry may be needed to forbid wca k at 
one loop [6a[28l, Chang, Keung (2004), Nelson (1984)]. 

This kind of Nelson-Barr type calculable 6> woa k can be 
mimicked in many extra dimensional models including 
superstring. Recently, for example # woa k is calculated to 
be 0(1O -12 ) from two loop level in a sequestered flavor 
and CP model [zl Cheung, Fitzpatrick, Randall (2008)]. 

Strictly speaking, the axion models also belong to the 
calculable models but we separate it from the models 
with spontaneous CP violation because there it is not 
needed to set 8 = 0. 



2. Massless up quark 

Suppose that we chiral-transform a quark as q —> 
e %1&a q. Then, the QCD Lagrangian changes as 

I d A x[~m q qq - 8{g 2 c GG}} -+ 

J d 4 x[-m q qe 2 ^ a q -(8- 2a){g 2 GG}] (11) 



where {GG} = (l/64 7 r 2 )e^'" T G° l/ G^. If m q = 0, it is 
equivalent to changing 8^8 — 2a. Thus, there exists a 
shift symmetry 8^8 — 2a. It is known that the tun- 
neling amplitude due to in stanton solutions with a zero 
mass quark vanishes J346I . 't Hooft (1976)], which im- 
plies that the shift symmetry is an exact symmetry. In 
this case, 8 is not physical, and hence there is no strong 
CP problem if the lightest quark (i. e. the up quark) 
is massless. The question for the massless up quark so- 
lution is, "Is the massless up quark phenomcnologically 
viable?" Weinberg's famous up/ down quark mass ratio, 
Z = m u /m d , gave Z = 5/9 J36J], Weinberg (1979)]. It is 
very similar to the recent compilation of the light quark 
masses, m u = 2.6^^ MeV,m d = 4.9±{\{ MeV, and 
Z = 0.48^^3 which is shown in Fig. [51 This compi- 
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FIG. 6: The allowed m u — nid region [268| . Two downward 
sloping lines are from the bound on [rn u + md)/2 and two ris- 
ing lines are from the bound on m u /rrid, which are determined 
by the masses of the meson octet. Two vertical and horizon- 
tal boundaries are from the Particle Data Book bounds on 
m u = [1.5, 3.3] MeV and m d = [3.5, 6.0] MeV 0. 



lation is convinci ng e nough to rule out the massless up 
quark possibility [2161 . Kaplan, Manohar (1986)]. In this 
review, we will use Z = 0.48 when a number is needed 
tho ugh t he appropriate bound may be 0.35 < Z < 0.60 
[53l. l268l . Buckley, Murayama (2007), Manohar, Sachra- 
jda (2008)]. 

For some time the massless up quark possibility was 
taken seriously [2161 . Kaplan, Manohar (1986)]. The 
reason is that even if the Lagrangian mass for the up 
quark is zero, the 't Hooft dctcrmincntal interaction may 
generate a useful up quark mass for the chiral pertur- 
bation. Ther e was a confusion on this issue for some 
time dSiH, Leutwyler (1990), Choi (1992)]. Now, it is 
clear that the massless up quark possibility is ruled out, 
even without using the l attic e calculation of the ratio, 
m u /m d = 0.410 ± 0.036 [Ml, Nelson, Fleming, Kilcup 
(2002)]. 



III. AXIONS 

The axion solution seems to be the most attractive 
one among three possible strong CV solutions, in partic- 
ular at present when the massless up quark possibility is 
excluded and calculable solutions need one loop suppres- 
sions. 



Peccei and Quinn tried to mimic the symmetry 9 — + 
9 — 2a of the massless quark case of (fTT|) , by con sidering 
the full electroweak theory Lagrangian J294I . 12951 . Peccei, 
Quinn (1977)]. They found such a symmetry if H u is 
coupled only to up-type quarks and H d couples only to 
down- type quarks, 

C = -q L u R H u - q L d R H d - V{H U , H d ) + (h.c.) - 9{GG}. 

(12) 
Certainly, if we assign the same global charge under the 



75 transformation to H u and H d , q 



e if} H u ,H d 

to 



*q, H u 



e l/3 H d , the flavor independent part changes 



C 



q L e-^ a u R e* p H, 



t/3 j 



9l6 



Y d R e if3 H d 



J/3 



- V{e^H u , e*H d ) + (h.c.) - (9 - 2a){GG}. (13) 



Since the full Lagrangian must possess the global sym- 
metry, the potential V should not allow the (H u H d ) 2 
term. Choosing j3 = a achieves the same kind of 9 shift 
of the massless quark case, which is called the PQ global 
symmetry U(1)pq. Unlike the massless up quark case, 
here 9 is physical. Even though the coefficient of {GG} 
changes in the same way in Eqs. (|lip and (fT3"|) , these two 
cases differ in that the tunneling amplitude vanishes with 
a massless quark (for which a detailed discussion will be 
presented in Subsec. IIII B[) and does not vanish without a 
massless quark. The reason is that the Higgs fields trans- 
form under U(1)pq, and one of the Higgs fields, called 
'axion' a, has the shift symmetry a — > a + (constant) 
and corresponds to the G oldst o ne b oson of the spon- 
taneously broken U(1) PQ [360|, [363, Weinberg (1978), 
Wilczek (1978)]. So we call the resulting axion from Eq. 
(fH)) the Peccei-Quinn- Weinberg- Wilczek (PQWW) ax- 
ion. If the consequence of the determinental interaction 
is only Fig. [TJ then out of two bosons rf and a, only 
rf obtains mass by the RHS diagram of Fig. [1] and a 
remains massless. If a remains mas sless , the strong CP 
problem is solved as envisioned in [294j, Peccei, Quinn 
(1977)] since for any 9 we can choose the VEV (a) such 
that the final 9 is zero. This was Peccei and Quinn's 
idea that (a) has a shift symmetry mimicking the mass- 
less quark case. However, a has interactions and can be 
produced in the stars and K meson decays, which differs 
from the massless quark case. 

At the classical Lagrangian level, there seems to be no 
strong CV problem. But, the axion coupling to {GG} is 
generated at one loop level, which is the U(1)pq-QCD- 
QCD anomaly. The 't Hooft determinental interaction 
we mentioned above is exactly this anomalous coupling. 
With this one loop term, the Lagrangian is not invariant 
under the phase shift symmetry /3, or a — » a+ (constant). 
Since it is explicitly broken at one loop level, the phase 
field P of the Higgs fields, or axion a does not have 
a flat potential, i.e. Fig. Q] is not complete. Wein- 
berg and Wilczek interpreted this phenomenon using the 
spontaneous symmetry breaking of the global symme- 
try U(1)pq. It is said that 9 is made dynamical where 



9 = a/F a , but in the PQWW axion case the component 
was there from the beginning in the phases of the Higgs 
doublet fields. The free energy depending on — cos 9 is 
the potential for the axion. Since it is proportional to 
— cos 9, the minimum o f the potential is at 9 = in 
CV conserving theories |353l . Vafa, Witten (1984)], and 
thus the vacuum chooses 9 = 0. We must discuss this 
effect below the chiral symmetry breaking scale as will 
be discussed in Subsec. If If Bl Thus, the axion solution 
of the strong CV problem is a kind of cosmological solu- 
tion. Note however that the weak C V v iolation shifts 9 a 
little bit, leading to 9 ~ 0(KT 17 ) [HI, Georgi, Randall 
(1986)]. 

The PQWW axion is ruled out quickly [H3, HHI, Don- 
nely et. al. (1978), Peccei (1978)], which was the reason 
for the popularity of calculable models in 1978 as dis- 
cussed in Subsubsec. Iff B II Nowadays, cosmologically 
considered axions are very light, which arises from the 
phase of SU(2)xU(l) singlet scalar field a. The simplest 
case is the Kim- S hifm an- Vainstcin- Z akharov (KSVZ) ax- 
ion model [Hi, [33l|, Kim (1979), Shifman, Vainstein, 
Zhakharov (1980)] which incorporates a heavy quark Q 
with the following coupling and the resulting chiral sym- 
metry 

C = - Q L Q R o + (h.c.) - V(\a\ 2 ) - 9{FF}, (14) 
C -> - Q L e l ^ a QB.e ll3 a + (h.c.) - V(\<r\ 2 ) 

-{9-2a){GG}. (15) 

Here, Higgs doublets are neutral under U(1)pq. By cou- 
pling a to H u and Hd, one can introduce a PQ symme- 
try also, not introducing heavy quarks necessarily, and 
the resulting axion is call ed the Di ne-Fischler-Srednicki- 
Zhitnitskii (DFSZ) axion [HiizSl, Dine, Fischlcr, Srcd- 
nicki (1981), Zhitnitskii (1980)]. In string models, most 
probably both heavy quarks and Higgs doublets con- 
tribute to the a field couplings. The VEV of a is much 
above the electroweak scale and the axion is a very light 
axion. 2 The SU(2)xU(l) singlet a field may mix with 
the Higgs doublet component by a tiny amount, whence 
practically we can consider the axion as the phase of a 
singlet field a, a = [(v + p) / \/2]e la ^ s with the identifica- 
tion a = a-\-2i:N£)wF a with the axion period 2itNE)wF a . 
Note that we use fs for the VEV of a or the value rel- 
evant in the field space and F a defined from the coef- 
ficient of the anomaly term. Namely, the coefficient of 
the anomaly {GG} defines F a as 9 = a/F a while the 
VEV(v) of er, a oc e la ' v , defines fs- The periodicity 2ir 
of 9 implies that F a cannot be larger than v = fs, and 
we have F a = fs/Now- It has been shown that mod- 
els with Now 7^ 1 has an ene rgy crisis problem in the 
standard Big Bang cosmology [3321 . Sikivie (1982)]. But 



models with Ndw — 1 do not have such a problem due 
to the mechanism to convert the two dimensional axionic 
domain wall disks surrounded by the axionic strings into 
radiation [U Barr, Choi, Kim (1987)]. 



A. Axion shift symmetry and reparametrization 
invariance 

In the original PQWW axion model, the Lagrangian 
in the effe ctive field theory language has been extensively 
discussed jl47l . 12931 Donnelly, Freedman, Lytel, Peccei, 
Schwartz (1978), Peccei (1989)]. Here, due to the sim- 
plicity in the formulae, we present the variant-type axion 
models where the PQ ch arg es are assigned only to the 
right-handed quark fields [25J, Bardeen, Peccei, Yanagida 
(1987)]. This discussion will make it easier in introduc- 
ing our general formulae below. The PQ current is [25l . 
Bardeen, Peccei, Yanagida (1987)], 

Ng N 

J ^ Q ^Fad^a + x ^2 dml^dm + (1/x) ^ ^Ril^ u m 

N g 

+ (-x) ^2 URi^um (16) 

i=N+l 

where N g is the number of families, N is the number 
of up-type quarks coupled to H u , and x = (H u )/(Hd). 
The color anomaly is non- vanishing, i.e. the divergence 
of JPQ is 



9 M JPQ = _ N{x 



1 a 

-^G^G^ 

X 'VK 



+m u u[i l5 e la ^' F « x ]u + m d d[i lb e la ~<'- x/F «]d (17) 

where we considered the one family model of u and d 
with N = 1. If TV were zero, there is no color anomaly. 
For a nonvanishing TV, we have to pick up the compo- 
nent orthogonal to the longitudinal Z^. Since the axial- 
vector part of Z^ current is proportional to J^ 3 , any ax- 
ial U(l) current orthogonal to the longitudinal Z^ is an 
SU(2)fl avor singlet current constructed in terms of right- 
handed quark fields. These include the currents corre- 
sponding to both rf and the PQ phase. Since 77' is known 
to be heavy, we integrate out 77' to obtain light fields 
below the chiral symmetry breaking scale. This corre- 
sponds to picking up an anomaly-free piece, orthogonal 
to the longitudinal Z^. It is 



ja _ tPQ 



l N ^ + bi 



1 



Z 



(m7m75« + Zd-y^d) 



(18) 



2 Once it was called an invisible axion |284| . |368| . Wise, Georgi, 
Glashow (1981), Nilles, Raby (1981)] but is better to be called a 
very light axion due to the possibility of its detection. 



where Z = m u /md. The divergence of (|18p is propor- 
tional to m u md which must be the case for the particle 
orthogonal to rj . 

Below we use the typical axion model (JT4")) because 
of its simplicity in assigning the PQ charges whenever 



an explicit example is needed. It has the following 
U(l)p Q charges, T, 



Field 



cr Ql Ql 



1 



where a is the phase field 
The current corresponds to 



In this example, the axial-vector current for U(1)pq is 
Jjl = Qlp.l5Q + vd^a where a is 
of a = (v/V2)e ia/v . 
the charge flow which satisfies the current conserva- 
tion equation if the symmetry is exact. But the 
axial-vector current is in general violated at one loop 
by the anomaly @, H Adler (1969), Bell, Jackiw 
(1969)], d^Jl = {N Q g 2 c /32Tr 2 )G« l ,G a ^, or d 2 a = 

{N Q g 2 /32n 2 v)G a ^G a ^ + (m Q /v)Qi l5 Q with the Q 
number Nq, which shows that the axion interaction with 
the SM fields is only the anomaly term (plus the anoma- 
lous coupling with the SM gauge fields). Here and in Eq. 
(fT7|) . we explicitly write the QCD coupling g 2 , but in the 
remainder of the paper we absorb the gauge coupling in 
the gauge fields except in the experimental section, Sec. 
IVl This axion is the one settling 9 at zero; thus one needs 
the axion-gluon-gluon anomalous coupling for which the 
color anomaly of J^ should exist. This kind of symmetry 
r is called the PQ symmetry. 

Axions are introduced as the Goldstone boson degree 
of a spontaneously broken glo bal U(l)pp symmetry in 
renormalizablc gauge models [2261 . |294 Peccei, Quinn 
(1977), Kim (1985)] and/or as a pseudoscalar degree in 
a more fundamental theory where the axion interaction 
arises as a nonrenormalizable anomalous interaction in 
the effective low energy theory. The most compelling 
nonrenormalizable interaction was observed in the com- 
pactification of 10-dimensional superstring models [37l[ 
Witten (1984)]. Below, we will treat that the axion is 
present as a dynamical degree at the electroweak scale 
whether it arises from the spontaneously broken PQ sym- 
metry or from a more fundamental theory with a non- 
renormalizable anomalous coupling, and focus on QCD 
interactions containing the axion degree, a = 9F a . Then, 
let us collectively write the most general form of its inter- 
actions: the c\ term is the derivative coupling respecting 
the PQ shift symmetry, the ci term is the phase in the 
quark mass matrix, and the C3 term is the anomalous 
coupling or the dctcrmincntal interaction £dct, 

£0 = Ifs^Od^e - ^G%G a ^ + {q L ilpq L + q R ilpq R ) 



+ c 1 {d li 9)vf'm - (q L m q R e lC2 ° + h.c 
+ c 3 ^G" G a ^ (or£ dct ) 



(19) 
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leptons.i 



where 9 = a/ fs with the axion decay constant fs up to 
the domain wall number (fs = NnwF a ), q is the fermion 
matrix composed of SU(3) C charge carrying fields. When 



the singlet scalar fields are easier to discuss, we use fs, 
and when the anomaly term is easier to discuss, we use 
F a . Ae P tons,e i s the axion interaction with leptons. Ci, ci, 
and C3 are pre-given coupling constants below the ax- 
ion scale fs with the mass parameter m defined to be 
real and positive below the electroweak scale. Then, the 
determinental interaction can be used instead of the C3 
term. 



£ dct = -2- 1 ic 3 6»(-l 



,N, 



-ic 3 e 



K 3N f~ 4 



F>et(q R q L )+h.c. (20) 



where we multiplied the overall interaction by 9 in the 
small 9 region and require the periodicity condition, 
c 3$ = c 3^ + 27r. The periodicity can be accomodated 
automatically if we replace — 2~ 1 ic^9 by 1, but then we 
must add a constant so that it vanishes at 9 = 0. The sign 
is ch osen such that the potential is a minimum at 9 = 
J35l Vafa, Witten (1984)]. With the fixed phases, the c 3 
term is given from the QCD vacuum structure (j4|) which 
does not have any dimensional coupling. But the instan- 
ton physics necessarily introduces the instanton sizes and 
hence a kind of QCD scale K for the interaction respect- 
ing the chiral transformation property for a flavor singlet 
operator Cdct- We will use either the anomaly term or 
£dct- The 9 dependence of the form ([20]) is — C39sm(cs9) 
which has the parity symmetry 9 — > —9. The Fourier 
expansion satisfying these is 

-2 -1 c 3 6» sin(c 3 (9) = -2 _1 [1 - cos(c 3 0)] + ^ a n cos(nc 3 9) 



n 2l a r , 



HO- 



where the Fourier coefficients satisfy, YTr, 

Neglecting the n > 2 terms, we will use just the cos(cs9) 

dependence. 

In the defining phase Eq. (fT9|) , the PQWW axion is 
given by C\ = 0, C2 7^ 0, and c 3 = 0, the KSVZ axion 
by c\ = 0, C2 = 0, and C3 7^ 0, the model-independent 
axion |37ll . Witten (1984)] in superstring models by 
C\ = 0, C2 = 0, and C3 7^ 0, and the DFSZ axion by 
C\ = 0, c-2 7^ 0, and C3 = 0. In general, axion models 
from high energy will have C2 7^ 0, and C3 7^ 0, and the 
shift symmetry allows c± =/= in a different basis. For 
simplicity, we discuss Eq. (fT^)) for one flavor QCD first. 
For Nf flavors, both Cj and 9 are defined from Nf x Nj 
matrices in addition to the anomalous coupling and hence 
the axion is included in Tr# which also contains the 77' 
meson part of QCD. For Nf flavors, Ci9 must be replaced 
by TiCi9. For the following illustrative discussion, we re- 
fer to one flavor QCD, but in Subsec. IIIIBl in the axion 
mass estimation we present the full Nf flavor QCD result 
with the chiral symmetry breaking taken into account. 

For the case of axion mass, ci, C2 and C3 terms may be 
relevant, but only the combination 02 + C3 appears. This 
Lagrangian has a shift symmetry a — * a+ (constant), 
which reparametrizes the couplings between ci,C2, and 
C3. Explicitly, the axion field dependent changes of the 



quark fields q^ 
c\ — ► c\ — a, c 2 - 



a iaa{x) 



q L and q R 



-iaa(x) 



q R give 



C2 — 2a, C3 — > C3 + 2a, and it must give 
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the same physics, i.e. |l69l |227| . Kim (1987), Gcorgi, Tomaras, Pais (1981)], 



ipi 



[a(x),A a Jx 



;ci,c 2 ,C3,to,A Q cd] = T 1PI [a(x),A^(x);ci -a,c 2 - 2a, c 3 + 2a,m, Aqcd]- 

I 



(21) 




a 



FIG. 7: The case with Now = 3 where three vacua are dis- 
tinguished. 



the 9 anomaly term which breaks the PQ symmetry. The 
global symmetry is not broken by the derivative term 
which therefore cannot contribute to the axion mass. 
From the reparametrization invariance (|2ip . the combi- 
nation c-2 + C3 is the correct combination for the axion 
mass, which is shown explicitly below. This derivation 
is included with a more complicated expression in the 
SUSY extension, but we show the c 2 + C3 dependence in 
this supcrgravity framework because it is the underlying 
symmetry in many axion models. Some of the following 
discussion is derived from [93|, Choi, Kim, Nilles (2007)]. 



The reparametrization symmetry dictates the non- 
derivative couplings satisfying C2+C3 = (constant), which 
is one of the reasons that we use 9 = #qfd + #qcd = 
#0 + $wcak as a physical parameter in axion models. Usu- 
ally, transferring all couplings of axion to the coefficient 
of GG, the axion decay constant F a and 9 are defined. 
Instead, if we use fs (defined to be the VEV of the sin- 
glet Higgs field a), there exists the coefficient C3 defined 
in Eq. (|19| . The triangle diagrams may give an integer 
times 9 and the instanton potential comes to the original 
value by a 9 shift of 2ir/(c2 + C3), with c-i + C3 = Ndw 
not necessarily 1 in the pseudoscalar field space. Thus , 
this integer is called the domain wall number Now [332L 
Sikivie (1982)] 



Now = |C2 



C3 



Trr(/colorod)^(/colored) (22) 



where the trace is taken over all heavy and light quarks 
and £ is the index of SU(3) C representation of colored 
fermions and the PQ charge is given for the left-handed 
chiral representations. The height of the potential is 
0(Aq CD ) of the nonabelian gauge interaction, which is 
shown in Fig. [7]with the domain wall number Now = 3: 
the bullet, the square and the triangle denote different 
vacua. Two important properties of axions in CP con- 
serving theories are: 

(i) periodic potential with the period 2irF a where F a 
is defined in (TIT))) with F a = fs/Now, an d 

(ii) the minima are at a = 0, 2nF a , 4:TrF a , 

This determines the cosine form of the potential. There 
exist the axion mixing with quark condensations as we 
will discuss it in more detail later. 

The derivative coupling, i.e. the c\ term, can never 
contribute to the PQ symmetry breaking effect, espe- 
cially to the axion mass. This axion gets its mass from 



Supersymmetrization 

We will discuss the reparametrization invariance with 
the SUSY generalization. In the Af — 1 SUSY mod- 
els with chiral fields z, there are the superpotential W(z) 
and the gauge kinetic function f(z) both of which are the 
holomorphic functions of z. The superpotential gives the 
C2 term and the gauge kinetic function gives the anomaly 
term C3. The PQ invariant Lagrangian, the Ci part, has 
the shift symmetry under the shift of the axion super- 
multiplct: A — > A + i(constant). This derivative cou- 
pling must appear from the D terms in SUSY models, 
i.e. through the Kahler potential. The real Kahler po- 
tential K (z, z*) must respect the PQ symmetry in the 
form of A + A, 

K = K [A + A} + {Z q [A + A]q iq2 + h.c.} (23) 

where the t? components of the fields are implied and 
'q's denotes quark supermultiplcts, 



q = tp q + iflipq 



(24) 



with the anticommuting variable &. Here, we used 1? for 
the anticommuting Grassmann number since 9 in this 
review is reserved for the axion 9 = a/F a . 



B. Axion mass 

The axion mass arises from the anomaly coupling 9GG. 
In this subsection, first we show that only the C2 and C3 
couplings are the relevant ones for the axion mass, and 
then we present the axion mass in the broken phase of 
the chiral symmetry. With SUSY, the discussion is a 
bit tricky, because the axion remains massless due to the 
massless gluino (as in the massless up quark case with 
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a spontaneously broken PQ symmetry). For the axion 
mass, therefore SUSY must be broken and here one has to 
look at all supergravity terms how they contribute to the 
axion mass. Nevertheless, we have the master formula 
pTj) for the axion which must be valid even when SUSY 
is broken. In this regard, SUSY is not special for the 
axion mass but only the anomaly consideration is the 
chief constraint. Thus, the following discussion applies 
even without SUSY, but let us discuss the axion mass 
in detail with the SUSY generalization to include the 
gluino effects and hence the C\ type derivative couplings 
to matter (quarks) and gauginos (gluinos). 

We have noted that there exists the famous anomaly 
coupling of the rj' meson which is the mechanism solving 



the old U(l) problem of QCD. In addition to rf , the axion 
a is introduced in the anomaly coupling and hence one 
must consider the mixing of r/ and axion [23|, l2a, 12231 
Baluni (1979), Bardeen, Tye (1978), Kim, Kim (2006)]. 

The C3 term is the anomaly coupling of the ax- 
ion, and we normalize the anomaly as the coefficient 
of e al3 ~ )S ei a E2f3ki 1 k2f3. With this normalization, from 
e a 0"/sd a A 13 d 1 A s leading to — c a ^ J ^^1^2^26, the C3 
term anomaly is defined with A3 = 1. 

It can be shown that, using the Kahler pote ntial (1231) . 
the kinetic energy terms of fermions contain [l!4J . 12831 



Cremmer, Ferrara, Girardello, van Proven (1983), Nillcs 
(1984)] 



J 



Y, z q fc# + -^7^ + \y q ,^fi6-^) + E (~ x ^ x - \ b l 



A 7 ^ A 



(25) 



where B^ and Y Qlti come from the auxilliary components of real Kq and Z q , respectively. In terms of the real parts 
TZ and y of Kq and Z (redefined from Z q of ([231) ). we obtain 
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(redefined from Z q of ([23])), we obt 
1 fdK 8K \ (8K 
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^ d A - ^d A 
d~A 



2\dA u < 



mr^ 



^d»A 

dA ' 



,{dlny\ , 

2 ' ST I d » a 



(26) 
(27) 



where 



G = -3 In ( — - ) +ln|IU| 2 , K = - e - Ko ' 3 , Z q = e z ', 



K = n + iK i = n, z q = y + n = 



y ■ 



(28) 



The C3 term is an anomaly term. In addition to the C3 
term, the c\ and c 2 couplings via loops of Fig. [8] will also 
generate the anomaly terms. The derivative coupling, if 
it has ever to contribute to the axion mass, should con- 
tribute to the axion mass by the anomaly through loops. 
In Fig. [SI the couplings for the triangle diagrams are 
represented in terms of c\ and c 2 . In supergravity mod- 
els, we consider B^ and Yq tft couplings which are nothing 
but c\. Consider a fermion with mass m. The derivative 
coupling through Fig. [5] has the anomaly co uplin g by 
picking up the coefficient of e 01 ^ 15 £\ a £2pki~ f k2f3 [l69l . See, 
for example, Gcorgi,Tomaras, Pais(1981)], 



A\= I dx\ 



where 



1— Xl 



dx- 



-4:f(x 1 ,X2;g,k 1 ,k 2 ) 
m 2 - f(x 1 ,x 2 ;q,k 1 ,k 2 ) 



(29) 



2;„2 



/ = (xi+x 2 )(l-xi-x 2 )q +2xi(l — xi — X2)q-ki+x l k 



(B ln Y Q4l ) 



(B m Y c 



c-xliTIs 




C17/.75 




k\,£\ k 2 ,£ 2 ki,£\ k 2 ,e 2 

(a) (b) 

FIG. 8: The Feynman diagrams for generating anomalous 
8GG couplings from c\ for a fermion with mass m. For C2, we 
replace ci^^-ys by c 2 m7 5 . 



Also, the quark mass term of Fig. [5] gives 

A 2 = dxi dx 2 — — — -. 

Jo Jo m -J{xi,xr,q, ki,k 2 ) 

From Eqs. ([2"9"[) and ([50)) . we construct 



I'M- 



-Ai +A 2 = 



/ dxi I 
Jo Jo 



2 dx 2 = 1. 



(30) 



(31) 
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When we calculate the axion mass in the real and pos- 
itive quark mass basis as usual, the anomaly coupling 
(a/F a ){GG} coupling (including the loop effect) is the 
sole source to the axion mass. In this basis, and also 
in any basis due to the reparametrization invariant com- 
bination C2 + C3, we do not have to discuss the deriva- 
tive couplings toward the axion mass. But, even though 



the derivative coupling generates the anomaly, it being 
derivative does not contribute to the axion mass. 

For one flavor QCD, we can check the above statement 
explicitly using Eqs. (|29I30I31|) . In the following two 
limiting cases, the integrals are easily computed as, using 
Eqs. ([29 ]) and ([30 ]) . 



J 



Case (i) : m < A QCD : r 1PJ = J-^T% 7 fc 2/3 f C3 + 2ci + 0(^] 



Case (ii) : m > Aqcd '■ Tip/ 



16tt 2 

r 



e Q/37 ^ l7 fc 2 /3 (c 3 + c 2 + O(-^r) 



(32) 
(33) 



Consider the quark mass term and the one flavor determi- 
nental interaction with the quark condensation, {(Il<IrI ~ 
Aq CD e"' If. Then, the potential takes the form 

V = m(q L q R )e lC20 + h.c. + (c 3 + CxAi + c 2 A 2 ) {g6} . 

(34) 
For the anomaly combination c 3 + c\A\ + c 2 .A 2 , the 
reparametrization invariancc Eq. (|2ip transforms c 3 + 
C1A1 + C2A2 to c 3 + 2a+ (ci - a)Ai + (c 2 - 2a)*4 2 = C3 + 
c\A\ + C2A2 where (|31[) is used, i.e. it is reparametriza- 
tion invariant. 



from which the axion mass is given by 



A 2 



IC2 + C3I 
fs 



Again the axion mass depends only on the combination 
c 2 +C3. Heavy quarks above the chiral symmetry breaking 
scale c, 6, and t give the c 2 term and vectorlikc heavy 
quarks above the electroweak scale give the c 3 term when 
we wrote Eq. (fill)) just below the electroweak scale. 



Axion mass with light quarks 



For Case (i), we consider the light quark below the 
scale Aq CD . Thus, we have 



= mv 3 cos ( j - c 2 9 j + A QCD cos f (c 3 + 2ci)0 + -j 



for which we choose c\ = 0. [If we kept c\, we must 
consider the kinetic mixing of a and 77'.] Integrating out 

the heavy 77' field as %■ — ~c 3 -f- from the Aq CD term 
which is the larger one, we obtain 



V = mv A cos (c 2 + c 3 ) — 
JS 



from which we obtain 



mA QCD 



[fg +C3J 
fs 



(35) 



Quarks u 1 d, and s belong to this category. 



For Case (ii), the heavy quark docs not condense and 
integrating out the heavy quark gives 

V = A QCD cos f (C3 + C 2 )^-J 



In the real world, there exist three light quarks whose 
masses are much smaller than the QCD scale Aqcd, and 
therefore the axion mass has the form anticipated in Eq. 
(|35|) . Even though there are two light quarks the axion 
mass dependence has the form ^Jm, because of F a ^> f n . 
This is because of the way in picking up the le ading term 
from the anomalous determinental interaction [223l Kim, 
Kim (2006)] as depicted in Fig. M 

In fact, this is simply obtained just by n oting that the 
instanton interaction is a U(l) singlet [227l Kim (1987)]. 
Suppose we integrate out quark fields; then the quark 
mass parameters appear in the effective interaction as 
shown in the first diagram of Fig. [9] In this vacuum in 
a theory with a massless quark, the tunneling amplitude 
vanishes so that the strength of the first diagram must be 
proportional to m q . With three quarks, we can generalize 
it as l/(r7i~ 1 +rn~ [ 1 +m~ 1 ). Suppose that there are only 
gluons and a very light axion a at low energy. Integrating 
out heavy fields, we are left with the flavor independent 
coupling aGG. Here we are not considering 7/ even be- 
low the quark condensation scale. If quarks are added, 
the flavor singlet coupling aGG can be split to quark 
mass terms with a u oc x/m u , aa oc x/rrid, a s oc x/m s , 
etc., as if quarks are not integrated over, m u ULUne. la + 
mdd[ J dRe lad0 + ■ ■ ■ , which shows that the flavor sin- 
glet coupling is of order 0(a/F a ). Then, even below 
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m s A 2 

A 


/ 


W e '(<=2+c 3 )0 


4- 

m u A 2 


y \A 

m d A 2 




i(c2+c 3 )6 



i(c 2 +c 3 ) 






m„A 2 



m d A 2 



771.A 




i{c^+c 3 )Q 



+ 0(m 2 A 4 v 3 ) 



FIG. 9: The 't Hooft determinental interaction. • denotes the 
quark condensation and x denotes the insertion of the current 
quark mass. The diagram highlighted with yellow predomi- 
nantly contributes to 77' mass, and 0(m u rrid) is neglected. 



the chiral symmetry breaking scale, we have the PQ 
charges proportional to l/m q . With this definition of 
quark charges, the axion mass comes from integrating 
out GG, and is proportional to a u + ad + a s which is 
~ m u mdm s / [m s (m u + rrid) + m u md], which was shown 
first for the PQWW axion 0, Baluni (1979)]. This is 
true even in the heavy quark KSVZ type axion models. 
Even if the light quarks do not have t he same PQ charge 
as in some variant axion models [25|, Il95l . 12381 . 12491 . I292L 
Krauss, Wilczek (1986), Peccei, Wu, Yanagida (1986), 
Bardcen, Peccei, Yanagida (1987), Kim, Lee (1989), 
Hindmarsh, Moulatsiotis (1997)], the axion mass has the 
same final form due to the reparametrization invariance, 
which will be shown below. So, the axion mass formula 
we write below is quite general. 



However, there was an assumption in this statement: 
rf was integrated out. So, it is necessary to include rf 
to obtain a more accurate axion mass. The light mesons 
and axion interactions must appear from those of Fig. [5] 
In this framework, however, the flavor singlet condition 
must be invoked as a constraint. [This flavor singlet con- 
dition is the anomaly matching without 77'.] Along the 
way, we would like to see how the \fm dependence arises 
below the chiral symmetry breaking scale in the KSVZ 
model. 

In the presence of vector like heavy quarks, the heavy 
fields are integrated out, whose sole effect is encoded in 
the low energy effective theory as nonrenormalizable cou- 
plings suppressed by F a , e.g. in the anomalous C3 type 
couplings with the SM gauge bosons. It is assumed that 
the heavy quark does not condense, since the QCD cou- 
pling is very small due to the asymptotic freedom at the 
heavy-quark Compton wave length scale and there does 
not exist a strong force to attract heavy quarks. Be- 
low the heavy quark scale, there is no massless mesons 
composed of heavy quarks. Therefore, the general form 
of the axion interaction Eq. (|19[) is valid at low energy. 
Firstly, the determinental interaction has the same chiral 
symmetry behavior as that of the anomaly term, and the 
anomaly term is removed in favor of the determinental 
interaction to include rf explicitly Second, let us choose 
the basis where the u and d quark masses are real. Since 
the strange quark mass is known to be below the QCD 
scale, we must include the strange quark with real and 
positive mass also in the instanton interaction. For sim- 
plicity, 7T° and rf , arising from quark conden satio ns uu 
and dd with decay constants f n and f V '{~ /*■) |164| . Gell- 
Mann, Oakes, Rcnncr (1968)], arc considered explicitly 
but with the r\ meson frozen. The effects of heavy quarks 
are included in the C3 term. If we keep c\ , the kinetic mix- 
ing of mesons and axion are present, due to the PCAC 
relation (0| J^ i (x)\xaesoiij(k)) = —ik^ffe~ lk ' x 8 i: j where 
J<ii ~ Ql^lhTiq. This would modify the axion mass, and 
hence it is easiest to calculate the axion mass by choosing 
the reparametrization parameter a such that c\ = 0. In 
this basis, denoting 7r°, rf and a in terms of dimcnsionlcss 
fields, 6^ = n //^, 8,y = rf/f v > f 6 = a/F a , we obtain the 
following effective interaction below the chiral symmetry 
breaking scale, 



J 



£ = -mviuLurieW'+W+W - m d {d L d R y^ +e ^ + ^ + h.c. + £ dot (36) 

where £det is given in Eq. (l20l) 

£dct = (-l) N, R'- 5 Uu L iLR){d L d R ){s L s R )e t[2e ^- C3(,) + ■■■ + flavor singlet constraint") + h.c, (37) 



where K has the mass dimension arising from QCD the anomaly (|3l?|) with c\ = 0. Note that t he log, det 
instanton physics. The above form is consistent with form in the effective Lagrangian was used in [1441 . Il45l 
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13571 . 13691 . I370I . Veneziano (1979), Di Vecchia, Vencziano 
(1980), Witten (1979, 1980), Di Vecchia, Nocodemi, Pet- 
torino, Veneziano (1981)] from the 1/N C expansion con- 
sideration, but we use Q37p because of its simplicity in 
the diagramatic expansion. The sign of the first dia- 
gram inside the box in Fig. [9] is de term ined to be mi- 
nus without the weak CV violation [353l . Vafa, Witten 
(1984)]. The QCD vacuum with the flavor indepen- 
dence of light quarks without the determinental inter- 
action chooses m q (qq) = — \m g \v 3 and we choose the sign 
of all quark masses be posi tive so that (qq) = (uu) — 
{dd) = -v 3 [HO, [HI SH SH, Dashen (1971), Lan- 
gacker, Pagels (1973, 1979), Gasser, Leutwyler (1982)]. 
Eq. (|37|) is the instanton interaction of Fig. [9j which 



gives A , m u A 3 , md A 3 , • • • by many ways of closing quark 
lines, shown in Fig. [51 but here one must invoke the fla- 



vor singlet constraint. The dominant term is the second 
diagram highlighted as yellow, which is flavor singlet and 
is the main source for the rf mass. 

Now let us restrict to the two flavor case. For the ax- 
ion, the key diagrams are those in the second line of Fig. 
|U If there are more than one QCD axion, then 0(m u md) 
diagram will be important at the next level axion mass. 
Integration over the instanton size includes very large 
instantons which covers the chiral symmetry breaking 
range where mesons appear as dynamical degrees, where 
we will invoke the flavor singlet constraint. The effec- 
tive interaction Hamiltonian of 9 n , 9^ and 9 = a/fs can 
be written, using the reparam etrization invariance (|21j) 
with Nf = 3 and 77 fixed, as [HI [HI Huang (1993), 
Kim, Kim (2006)], 



- V — m u v 3 cos (6 V + 9 V >) + m d v 3 cos (-9„ + d v >) + -r^ cos (29 v > - (eg + 4 + c 3 )6) 



2„6 



+m l 



■ K 5 



(c« + c d 2 +c 3 )9)+m d ^- 



. (9„ + 9 V , - (c 2 " + 4 + c 3 )9) , 



(38) 



r 



where A u and A d are parameters describing the result 
of the Feynman and instanton size integrations. The 
(— 1)^ is cancelled by the fermion loop or (—v) factors. 
If m u = ma, A u and Ad are equal. For m u ^ md, A u 
and Arf must be different. The instanton interaction is 
flavor independent, which should be respected in the in- 
teraction (|38p . The m u and md linear terms from the de- 
terminental interaction should be flavor independent, i.e. 
m u A 2 u + rndAd = flavor independent. Since it vanishes if 
one quark is massless, it must be a function of m u md- 



Thus, the instanton size integration with current quark 
masses must give m u A^ + mdA d = 2m u mdL 2 / (m u + md), 
which vanishes if any quark is massless. This is because 
the original gluon anomaly term {GG} does not distin- 
guish flavors, and the smallncss of the current quark 
masses enable us to expand the 't Hooft determinental in- 
teraction in terms of powers of the current quark masses. 
Then, the 3x3 mass matrix M 2 of a, 77' and 7r°, tak- 
ing into account of the chiral symmetry breaking and the 
solution of the U(l) problem is given as 



M. 



a,r)' 



1 

( c 2 [A?,, + 2nAl st }/F 2 -2c[A*, + fi^J/f'F 

-2c[A*, + nAl st ]/f'F [4A* + 2MLt + m + v 3 ]/f 2 -m.v 3 /ff 

-m-v 3 Uf (m + v 3 + 2 M A 3 nst )// 2 J 



(39) 



where c = c£ + 4 + c 3 , F = fs, f = U, f = U, A*, = limit f/F, f'/F « 1, we obtain 



v 6 /K /2 , and Af nst = L 2 v 3 /K 2 ,m+ 
md — m u , and 



m u + m dl m^ 



M = 



m u m d 
(m u + m d ) 



(40) 



Certainly, Eq. (J35]) realizes the solution of the U(l) prob- 
lem due to the A*/ term in the (22) component. In the 



2 m + v 



- 2ML 



f 2 

■/ 7T 



2 4Af , + m+v 3 + 2f,Af nst 

m n , ~ '- 



■hi' 



ml ~ 



F 2 (1 + Z) 



;f 2 mlo (1 + A) 



(41) 
(42) 
(43) 
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where 






(44) 



In this form, the 7r mass has the standard m+t> plu s the 
instanton contribution to light quark masses [871 . 12161 
Kaplan, Manohar (1986), Choi, Kim, Sze (1988)]. From 
(|41l42p . we estimate the parameter Ai which is the 
source of the solution of the U(l) problem: A*, = 

ifv' m v' ~ /' m ')/ 4 ~ ( 202 McV ) 4 with U ~ 86 McV 
and fn ~ 93 MeV. In any axion model, this form is valid 

with \c\ = Nr>w- Using the standard definition on the 

axion decay constant F a = F/c, we obtain 



f>l 



d + zy fi 



;i + a). 



(45) 



So, even though the instanton diagrams of Fig. [9] con- 
tain the linear quark mass diagrams summed over, the 
diagonalization process with mesons signals the predom- 
inant contribution of the lightest quark. The flavor sin- 
glet condition we discussed before chooses the following 
linear quark mass dependence 



/; 



1 
m,. 



1 
m d 



(46) 



Neglecting instanton contribution to the current quark 
masses, we obtain m a ss 0.60 cV [10 7 GeVjF a ], for the 
mass ratio Z ~ 0.48 as summarized in 268, Manohar, 
Sachrajda (2008)]. An earlier frequently cited Z is 5/9 
dUGM Weinberg (1979), Gasser, Leutwyler (1982)]. 
The correct axion mass has to include the current quark 
mass change due to instantons. However, the resulting 
estimate of A turns out to be small. 



Comparison with old calculation 

Now, let us comment on the old anomaly matching 
condition. If any quark mass is zero, there exists an ex- 
act symmetry a — > a + (constant), i.e. axion is massless, 
above the chiral symmetry breaking scale. Below the 
chiral symmetry breaking scale, it is likely that this con- 
dition is satisfied. Let us denote the original current as 
Jpq . This current is anomalous above the chiral symme- 
try breaking scale, <9 M J£ Q = (N Q /S2Tr 2 )G a ^G a ' lv where 
Nq is the number of heavy quarks with T = 1/2. Below 
the chiral symmetry breaking scale, we considered two 
pseudoscalar mesons which have anomalous couplings: 
rf and a. The global anomaly matching conditi on w ill 
work if there were no chiral symmetry breaking [3471 't 
Hooft (1979)]. For chiral symmetry breaking, there is 
no massless fermions and we consider only color singlet 
mesons below the chiral symmetry breaking scale. Thus, 
the bosonic current must be anomaly free after integrat- 
ing out all heavy fields including rf , i.e. we consider an 



anomaly free current J£ instead of JpV below the chiral 
symmetry breaking scale J2271 Kim (1987)], 



t/j, _ jf* _ 



Nr 



2(1 + Z) 



(w7 M 7 5 u + ZdriyTihd) (47) 



where the divergence of the second current gives a singlet 
pseudoscalar density so that axion does not mix with ir°. 
The expression (|45p with A shows that the finite rf mass 
enters into the a — rf mixing. 



Mesons without axion 

Even if there is no axion, we can diagonalize the mass 
matrix. If m u = 0, one starts with an exact up quark 
chiral transformation, which would lead to a Goldstonc 
boson 9 in the vacuum, (uu) ^ 0. This Goldstone boson 
couples to neutron through (ci 9^0)777^7577. In reality, rf 
obtains mass by the anomaly, and the symmetry remains 
unbroken: it is the phase symmetry of (uu). Therefore, 
any violation of the shift symmetry must be invoked such 
that it g oes a way in the limit fj, — ► 0, which is Dashen's 
theorem [HI Dashen (1971)]. Thus, from Eq. fl35J) we 
obtain the VEVs of rf and ir° for a small 9 



vV) 



/ 



a , . uv° 

-2 {1 + Z) k 



171+ 



(48) 



The VEVs of rf and ir° are vanishing if 6 — or any quark 
mass is zero. In addition, we can estimate the rf proper- 
ties from the interaction (v 9 / K 5 ) cos(2rj' / f v ') where f v i 
is the ?y' decay constant and K has a mass dimension. 
It comes from the yellowed diagram of Fig. [9] Compar- 
ing 7T° — > 27 and 77' — > 27 decay widths, 7.74 eV and 
4.3 keV, respectively [9|, Amsler et. al. (Particle Data 
Group, 2008)], we obtain /*, = (4/3)(m^/m^.)(r(7r -> 
2 7 /r(77' -> 2 7 ))/2, or ,/V « 86 MeV. Fitting to the rf 
mass, we obtain K = (v 9 / f*,m'*) 1/5 = 240 MeV. 



The 8 — vacuum with axion 

We have shown above that the Lagrangian ([38]) chooses 
9 = in CV conserving theories if 9^ = and 9„i = 0, 
which is determined by QCD dynamics. However, if 
CV is broken, the vacuum value of 9 is shifted from 
the 9 = value, by the presence of any linear term of 
9 V ,9 V ' or/and 9. The meson potential is invariant un- 
der CV with CV{n) = CP{rf) = CV{a) = -1. Such 
linear terms are generated by considering CV violating 
phases and chirality flipping (L<->R) insertions. So, lin- 
ear terms of a are generated by combining the 't Hooft 
determinental interaction and CV violating weak inter- 
actions. Linear terms of 7r and 77' can also be generated 
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by considering the weak interactions alone without the 
determinental interaction, but the conditions of the fla- 
vor singlet, chirality flipping (L<->R) and CP violating 
effects do not occur at one loop level. In th e SM with 
the Kobayashi-Maskawa CV violation, Ref. [1531 . Ellis, 
Gaillard (1979)] shows that a finite correction occurs at 
the fourth order, 0(a 2 ), leading to a very small NEDM, 
but infinite corrections occur from 0(a 7 ). These can 
give rise to a linear term of n. In the SM, the pioneering 
calculation with axion has been perfor med in ch iral per- 
turbation theory to obtain < 1CT 17 fl67L EH, Georgi, 
Kaplan, Randall (1986), Georgi, Randall (1986)]. The es- 
timated however is far below the current experimental 
limit of 1CT 11 . 



In the barred notation, there exist only c\ and c 2 . 

We will discuss one family without separating ci ]2 
into c"' 2 first for an illustration, and then we will dis- 
cuss the cases with c"' 2 and write down formulae for 
three families. Let us define the initial parameters c\ , c 2 
and C3 together with the definition of the vacuum an- 
gle 9q = Oqcd- In principle, the initial vacuum angle 
can be a free parameter. Here, the vacuum angle Oqcd 
is defined such that c\ = 0. Picking up the axion de- 
pendent chiral rotation charge defined below the chiral 
symmetry breaking scale Eq. (|47[) . the chiral quarks 
in the chiral perturbation theory are transforming as 
q L -» cxp{iQ A &)qL,qR -> exp(-iQ A 0)q R where 



C. Axion couplings 

The axion interactions arc given in Eq. (|19p which 
are depicted in Fig. ITOl where we have not drawn aWW 
and aZZ diagrams which are orthogonal to the 077. The 
diagrams of Fig. [10] arc complete for the low energy axion 
phenomenology, where the suppression factor, 1/F a , by 
the axion decay constant is explicitly shown. 




FIG. 10: The Feynman diagrams of axion couplings. G and 
7 are gluon and photon, respectively. C3 and c a77 couplings 
are anomalous. 



Qa = 



1 M- 



2 TrM" 1 ' 



M' 1 =diag.( 



1 



J_. 



(50) 



The derivative interactions of axion is obtained in this 
way |l67l . 12151 . Kaplan (1985), Georgi, Kaplan, Randall 
(1986)]. 

For the KSVZ axion, we have c\ = c 2 — and C3 = 1, 
and the coefficient of the gluon anomaly term is (a/F a ) + 
Oqcd- Hence, redefining the axion as a + F^qcd, we 
obtain 3 



KSVZ axion (c x = 0, c 2 = 0) 

Cl = 2-C 3 = i, 
c 2 = c 2 +c 3 = 1. 



(51) 



Here, c 2 must be split according to the flavor singlet con- 



dition to c 2 



For the DFSZ and PQWW axions, a = 0, c 2 ^ and 
C3 = 0. If a non-vanishing Oqcd is introduced here, we 
have, using the rcparametrization invariance (j2 1[) . c' x — 



— c 2 /2,c 2 = 0, and c 3 = c 2 . 

the gluon anomaly term is 02(0,/ fs) 

redefining axion as a + (/s/c 2 )6*qcd and going back to 

the C3 = basis, we obtain for one family 



Then, the coefficient of 
qcd, and hence 



1. Axion hadron coupling 

When we discuss axion-hadron interactions, which are 
the relevant ones for low energy laboratory experiments 
and physics at the core of supcrnovae, we must integrate 
out gluon fields. Technically, it is achieved by using the 
reparametrization invariance to remove the C3OGG cou- 
pling. If we keep the C3 coupling, we must consider the 
axion-gluon-gluon interactions also, which is hard to be 
treated accurately with its face value but must be the 
same as the one in the C3 = basis. In this way, the 
quark interactions are changed from the original value as 



(49) 



ci - 


-> Cl 


= Cl + \c z 


C2 - 


■+ c 2 


= C2 + C3 


C3 - 


-* c 3 


= c 3 - c 3 = 



DFSZ and PQWW axions 

Cl = \{~c 2 +c 2 ), 

c 2 /o, ca = 0. 



(52) 



Again, c 3 must be split according to the flavor singlet 
condition to c 2 + c 2 according to the anomaly matching 
condition, Eq. ((47)) . 

Integrating out the heavy a field and heavy quark 
fields, the massless (at this level) degree a = F a which 
appears from the phase of the singlet field a = ((a) + 
-%)e ie appears in the effective low energy Lagrangian. 
If there are multiple SM singlets Si carrying PQ charges 



The sign convention is stated below. 
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and VEVs, then the axion component is proportional to 

o = ^£ I WS0' K = (]Tr^ 2 ) 1/2 (53) 



V a 



where a(Si) is the phase field of Si. The PQ charges are 
defined such that the smallest nonzero absolute valuc(s) 
of the PQ charges is 1 so that every scalar field returns 
to its original value after 2tt shift of its phase. Let us 
discuss axion couplings after integrating out the heavy 
fields. 



In the KSVZ model C3 is calculated by the 
triangle diagram of heavy quarks for the global 
anomaly. The domain wall number Now is Now = 
Trr(Q L );(Q L ), with F a = V a /N DW where T(Q L ) (de- 
fined as Q L -> e ir( - Q ^ e Q L under a -> a + F a 6) is the PQ 
charge and 1(Ql) is the index of SU(3) C representation. 
Every field is represented in terms of left-handed fields, 
and the PQ charges are defined such that the SM sin- 
glet a coupling to heavy quarks carries one unit of the 
PQ charge. If the light quarks also carry the PQ charge, 
then ([22)1 gives Now-, which belongs to the generic very 
light axion model discussed below. The anomaly calcu- 
lation gives the one loop coupling (N ow a /Va){GG} : but 
since the vacuum angle 8 or axion is given by the coef- 
ficient of {GG}, F a is defined by dividing V a of (|55|) by 
Now and hence C3 = ±1 where the sign coincides with 
that of Ti:T(Ql)1(Ql). As a convention, let us choose it 
to be C3 = +1, which is choosing the effective PQ charges 
of heavy quarks to be positive. Transferring C3 to C2, we 
split C3 = C2 + c 2 using the PQ charges of l|50p. 



KSVZ axion 
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(54) 
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flavor singlet condition, we obtain for one family 



DFSZ axion for one family : 
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where v u = \(>/2H°)\,v d = \(V2H°)\,v EW 
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The PQ charges c^ 
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and 



W 



(55) 



(v 2 u + 
oiH u 



and Hd are obtained by considering the orthogonal com- 
ponent to the longitudinal mode of Z boson. Remember 



u.d 



that the signs of c 2 ' are chosen from the convention that 
the PQ charges of H u ,d are positive. This result is for one 
family. 



If we have N g families, we can calculate the couplings 
just below the electroweak scale where all quarks ob- 
tain masses. Then, we obtain for three families, c 2 = 



l«<*r 



and 



\Vu\" 



Us- 



ing the reparametrization invariance, we can calculate 



-1' L 2 



and c 3 , just above 1 GeV: 



0, cV = 



2 C 2 



and 



c 3 = N g (J2i c 2) ~ N g . Then, we integrate out heavy 
quarks c, 6, and t to obtain the effective couplings just 
above 1 GeV, which does not introduce any new C2 terms. 
Now, there are three light quarks u, d, and s for which 
we use the reparametrization invariance to remove the 
C3 term such that the isosinglct condition is satisfied, 
d^J? is anomaly free where J° = J^ Q — a. u u-f^-f5U — 

add^^d - 
a u + a d + a, 
condition of d^Jf, determines 

H 1 



a s sj^ 5 s and d^J^ = N g {GG}. Thus, 
1 = N g is satisfied and the SU(3)fl av0 r singlet 
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with 



In the DFSZ model, c 2 and c 2 are calculated by trans- 
ferring the phase of a to H u and Hd by the PQ sym- 



metry such that {H u 



^/2 V2e *r u a/v„ anc i 



m 



J2v x e lVdalV ° if H*H*cr 2 
in terms of PQ charges T u and Td of H u and Hd- 
Here, a = V a 6 is not the mass eigenstate and instead 
of V„ the mass eigenstate a uses the decay constant 
F a = \(T U + TdfV* + Tlvl + I>2]i/2 „ {Tu + Td)v<j 
for V a 3> v u ,Vd, and the axion component a = [(T u + 
^d)V a a + T u v E wa(H u ) + T d VEwa(Hd)]/F a ~ a, and 
F a = VajNow ■ So, in the DFSZ model they are given by 



if H^H^cr defines the PQ charge of a 



|5& Carena, Peccei (1989)]: 



_ WV 



\vuV 



and 



c 3 = 0. By the reparametrization invariance, Eq. 



(p?Tj) . we can use d{ 
and c' 
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( 2 
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lu 
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0, 



1. Removing c 3 according to the 



N q m u m d m s 



m u md + m u m s + mdm s 



N g (m u + m d )Z 

{i + zf ■ 



Therefore, removing c 3 by the reparametrization invari- 
ance, we obtain 



DFSZ axion for N g families : 
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J EW 



J EW 



(57) 
(58) 



If heavy quarks and also H u d carry PQ charges, we 
must consider all these. If one SM singlet a houses the 
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axion, then wc obtain 



General very light axion : 
;(1±JV») 
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(59) 
(60) 
(61) 
(62) 



where the PQ charges (+ or -) of H u and Hd deter- 
mine the sign (- or +) in front of DFSZ component and 
5h = 1 or if the corresponding Higgs doublets carry the 
PQ charges or not. For the MI axion from supcrstring 
which is a hadronic axion, in principle there can exist an 
additional contribution to c\ as will be pointed out in 
IVIF11 

If there is no heavy degrees carrying the PQ charges 
above the electroweak scale, then c^ in the so-called 



PQWW model is given by the PQ charges of H u and 
H r :, 



PQWW axion : Same as Eqs. (|57I58|I . 



(63) 



All models have Cj and c^. For the original c 2 term, 
different models give diff erent values, fo r example some 
variant axion models |25l . Il95l 12381 I249J Krauss, Wilczek 
(1986), Bardcen, Peccei, Yanagida (1987), Kim, Lee 
(1989), Hindmarsh (1997)] have different C2S from those 
of the PQWW axion. For astrophysical application, we 
must keep both c\ and C2- The c\ and C2 terms respec- 
tively give the axial vector and pseudoscalar couplings. 
The axion operator in the flavor SU(3) space can be writ- 
ten as 



K% " ci 2 )F 3 
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(64) 



where F$ and Fs/y/3 are the third component of the 
isospin and the hypercharge operators, respectively, and 
1 is the identity operator. 



J 



Then, the derivative couplings with nucleons and mesons below the chiral symmetry breaking are defined as 
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Here the axial vector coupling parameters of the nucleon octet are given by F = 0.47, D — 0.81, S ~ 0.13 ± 0.2 [9|, 
Amsle r et a l. (Particle Data Group, 2008)]. For example, for the hadronic axion couplings we obtain the results given 
in [nlHH, Kaplan (1985), Chang, Choi (1994)], 
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For the DFSZ axion, there exist additional contributions 
from the extra terms in Eqs. (|57I58[) . 

Similar expressions might be attempted for the pseu- 
doscalar couplings in terms of c^' and the pseudoscalar 
coefficients F',D' and S". But, for the axion current, 
corresponding to J" there does not exist an anomaly as 
discussed in Eq. (|47[) and we do not write down the axion 
pseudoscalar couplings. The anomaly carried by axion 
above the chiral symmetry breaking scale is left over to 



rf below the chiral symmetry breaking scale and hence 
these pseudoscalar couplings are for the r)' meson. The 
axial vector current of 77' to the nucleon octet N = q®q®q 
is, 



(70) 



JJ1 = fadyrf + g 5 N N^ l5 T N 

where To is properly normalized TrTg = 

l/y/2Nf, and g% is determined by strong interaction 



I or T 
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dynamics. The original global symmetry breaking term 
(j2"0")) is transferred to N '^N ne laiV ^v' which is actually 
the nucleoli mass term, 



AC = -m N N L N R e ia '^' /f i' + h.c. 



(71) 



For example, the SU(6) wave function of spin-up neutron 

is 






2d^u^S - 2SuU^ - 2uU i (f 

- 2<tfSvl - 2d i d}vl + 4d T dM) 

, (72) 
where the quarks are now interpreted as constituent 
quarks below the chiral symmetry breaking. At low en- 
ergy, the symmetry is the only relevant one for consid- 
eration. The octet charge a[ is determined by strong 
interaction dynamics. The ducaplet has a different U(l) 
charge a'{. Two anomaly matching conditions, PQ-B-B 
and PQ-Qcm-Qcm, may be used but do not give a useful 
information because of many form factors. So. simply 
interpreting the PQ charges of the current quarks being 
transferred to the constituent quarks in the octet with a 
multiplcation factor g b N , we obtain the PQ charge of neu- 
tron as the PQ charge of one constituent quark. Thus, 
NlNr has the phase a' = 2g b NX /2Nf/Nf. Guessing that 
g b N is similar to the octet form factor g^ ~ 0.75 [l66l 
Georgi (1984), p.100], a[ is estimated as 1.22. 



for the 20% allowance fr om t he tree level chiral pertur- 
bation theory estimation [216l Kaplan, Manohar (19 86)] . 
For illustration, let us take c x sb ~ 1-98 for Z ~ 0.5 |268t 
Manohar, Sachrajda (2008)]. 

In the KSVZ model, c a is determined by the PQ 
charge carrying heavy fcrmions. If there is only one neu- 
tral quark for this, then c° would be zero. If there 
is only one PQ charge carrying heavy quark with the 



electromagnetic charge Q em , then c a 



Q. 



2 

em* 



But, 



realistic models from a fundamental theory it is more 
likely that there exist many PQ charge carrying quarks, 
and the coupling given for one PQ charge carrying heavy 
quark is presented just an illustration. 

In the DFSZ model, we consider only light quarks and 
leptons. The PQ charges of H u and Hd determine the 
PQ charges of u and d quarks. For the PQ charge of e, 
we have two possibilities: Hd gives mass to e and the PQ 
charge of e is the same as that of d, or H u gives mass to 
e and the PQ charge of e is the opposite to that of u: 
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(78) 



2. Axion-photon-photon coupling 

As we calculated the c$ coupling for the KSVZ axion, 
we can calculate the axion-photon-photon coupling, by 
simply substituting the gluon lines by photon lines and 
the quark triangles by charged fermion triangles. Since 
we are interested in low temperature experiments, we 
consider the energy scale below the electron mass. There- 
fore, considering V a = Nr>wF c 
PQ charges of charged fermions becomes 



c" 77 calculated from the 



Trr(Q L )Q 
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em 
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Now 



(73) 



Below the QCD chiral symmetry breaking scale, we chiral 
transform light quarks to obtain 



C u 






(74) 



where 



where the PQ charges of H u ,d were chosen to be positive 
before. So, in applying Eq. (|75p . we must choose the 
PQ charges of light quarks to be positive and hence the 
signs of (J77I78P must be reversed. For the PQWW axion, 
the coupling is the same as those of (J77I78P with positive 
signs. 

For a general light axion, the axion-photon-photon 
coupling depends on the ultraviolet completion of the 
theory If the axion mass is lighter than 2m e , its lifetime 
is 
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(79) 



where we used Z ~ 0.5 and the age of universe tjj s=s 
4.35 x 10 17 s. For c a77 = 0(1) the axion with 24 eV 
mass has the lifetime tjj [l90t 12761 . Moroi, Murayama 
(1998), Hannestad, Mirizzi, Raffelt, Wong (2008)]. 
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where the chiral symmetry breaking effect is, including 
the strange quark mass effect, 



-xSB 



[(4 + 1.05Z) 
1 + 1.05Z 



[1.762,2.260] 



(76) 



3. Axion lepton couplings 

The tree level axion lepton (l) coupling arises in the 
DFSZ and PQWW axions where the lepton mass term 
through the PQ charges of Hd or H u defines the c; of Fig. 
1101 The removal of C3 term does not change the coupling 
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Ci , and hence we obtain the following tree level couplings 
of axion and lepton: 



DFSZ 
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NdwFhV ew 
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wF a v 2 EW 
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N 



lijsla, lepton mass by Hd (80) 



li'ysia, lepton mass by W u (81) 



where the PQ charges of H u .d arc chosen to be posi- 
tive. For the PQWW axion, just F a is replaced by vew- 
For the KSVZ axion, the axion lepton coupling occurs 
at higher order which is negligible in astrophysical appli- 
cations. For the generic very light axion, the couplings 
given in Eqs. (|80I81|) arc applicable. 

Even though, the tree level coupling of the axion with 
electron is absent in the KSVZ model, the axion-electron 
coup ling is present at one loop through the c a77 coupling 
pll Srednicki (1985)] 
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where A is the chiral symmetry breaking scale and N^ w 
must be multiplied in models with Njjw 7^ 1- On the 
other hand, the DFSZ axion coupling to electron is 
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where N g = Ndw/2 is the number of families and X^ 



sin 2 /3 



v u/ v ew f° r t ne case 0I Eq. ((80 



D. Old laboratory bounds on F a 



With the axion couplings discussed in Subsec. IIIIC1 
one can estimate the axion production rates in vari- 
ous experiments. The null experimental results give the 
bounds on the relevant axion couplings. These ar e ex- 
tensively discussed in earlier reviews [73, 12271 1293| . Kim 
(1987), Cheng (1988), Peccei (1989)]. These old labora- 
tory bounds, immediately studied after the proposal of 
the PQWW axion, basically rule out the PQWW axion, 
i.e. give the axion decay constant F a greater than 0(10 
TeV), 

F a > 10 4 GcV : old laboratorybound. (84) 



IV. AXIONS FROM OUTER SPACE 



From Eq. ([79"]) , we note that the axion lifetime is longer 
than t\j for m a ~ 24 eV, and this kind of axion is impor- 
tant in cosmology. For m a < 23 kcV with c a77 — 1, the 
axion lifetime is longer than 10 min, allowing solar gen- 
erated axions below this mass have time to reach Earth. 



These examples illustrated the importance of studying 
low mass axion effects in astrophysics and cosmology. 

The window for F a obtained from the astrophysical 
and cosmological constraints is given by 



0.5 x 10 9 GeV <F a < 2.5 x 10 12 GeV 



(85) 



where the upper bound is understood with the initial 
misalignment angle of order 1 . 



A. Axions from stars 

In this Subsec. we present the key arguments leading 
to the axion constraints from astrophysical sources. Ax- 
ions have very small masses and therefore can be emitted 
without important threshold effects from stars, in anal- 
ogy to neutrinos. The method to constrain axion models 
is basically the overall energy loss rate, whether using 
the individual stars (e.g. Sun, SN1987A) or the statisti- 
cal properties of stellar populations (e.g . the sta rs in a 
globular cluster being a test population) [246l 13061 . Kolb, 
Turner (1990), Raffelt (1996)]. 

We may use the axion couplings to j,p,n, and e to 
study the core evolution of a star. The simple bounds are 
obtained just by comparing the energy loss rates by the 
axion emission and by the neutrino emission. Studying 
the evolutionary history of a star by the axion emission 
may give a stronger bound than the one obtained from 
the energy loss rate but may not be so reliable compared 
to the one obtained from the energy loss ra te. Since there 
are good reviews on axion astrophysics [2, 1304 I307J . l35ll 
Turner (1990), Raffelt (1990, 2008), Amsler et al. (Parti- 
cle Data Group, 2008)], here we briefly comment on axion 
physics in stars (Sun, low mass red giants, supernovae) 
to cite the reliable F a bound. 

With the axion emission, the Sun consumes more fuel 
and needs an increased core temperature. From the Pri- 
makoff proces s 7 + Ze — > a + Ze in the hadronic axion 
models, Ref. ^M, Schlattl, Weiss, Raffelt (1999)] gives 
the axion emission rate L a ~ 3.7 x 1O~ 2 L with a 20% 
increase of the 8 B flux with the increased core temper- 
ature. The 8 B neutrino flux gives the best bound on 
the solar axion emission rate. The measured 8 B neu- 
trino flux 4.94 x 10 6 cm" 2 s- 1 @, i, SNO Collaboration 
(2005a, 2005b)] is consistent with the axion emission if 
L a < 0.04L Q (IH Bahcall, Serenelli, Basu (2005)]. This 
translates to the F a b ound of F a jc ail > 2.6 x 10 6 GeV 
for L a < 0.04L o [ill, Schlattl, Weiss, Raffelt (1999)]. 

For the axion-electron coupling such as in the DFSZ 
axion models, the axion emissio n fro m the globular clus- 
ters gives a useful F a bound J309I . Raffelt, Dearborn 
(1987)]. Stars in a globular cluster are assumed to have 
the identical Y (helium fraction) and metallicity frac- 
tion. The helium core before ignition is degenerate and 
the bremsstrahlung emission is very effective, whereas 
the Primakoff emission is suppressed by the large plasma 
frequency and the helium ignition docs not give a use- 
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ful F a bound for the KSVZ axion. However, after he- 
lium ignition the core degeneracy is lifted, the Primakoff 
effect becomes important, and the consumption of he- 
lium fuel is accelerated by the axion energy loss, short- 
ening the helium- burning lifetimes. Horizontal branch 
stars in several globular clusters confirm the expected 
helium-burning lifetimes, which agrees with the stan- 
dard prediction and the axion losses should not exceed 
e a < 10 e rg g _1 s _1 in the cores of horizontal branch stars 



[64j, |305j, Raffelt (1990), Catelan, de Freista Pacheco, 
Horvath (1996)], which leads to F a /c ail > 2 x 10 7 GeV, 
a factor of 10 improvement over the solar bound. Note 
that this globular cluster bound is for models with an 
appreciable axion-electron coupling. 

In the study of the axion emission in the small mass red 
giants, the processes j+Z — > a+Z, e + Z — > a+e + Z, and 
7 + e — » a + e were considered. The early studies were the 
simple co mparison of th e axion emission and the neutrino 
emission J159I . Il60l l250l Fukugita, Watamura, Yoshimura 
(1982a 1982b), Krauss, Moody, Wilczek (1984)]. In the 
study [124l Dearborn, Schramm, Stcigman], it is summa- 
rized as F a > 2.1 x 10 7 c a77 GeV if the Primakoff process 
7 + Z -> a + Z dominates and F a > 3.7 x 10 9 sin 2 /3 
GeV if the Compton process dominates (which is for the 
DFSZ axion, viz. Eq. (|80"j)). The Primakoff process is 
present in any axion model, and hence the Primakoff pro- 
cess bound is almost model-independent except in the 
region m a > 200 kcV where a is too heavy to be pro- 
duced in the core of a star. But this threshold effect 
is irrelevant since the PQWW axion region is excluded 
already. Note however that there is no confirmed obser- 
vation of neutrinos from the small mass red giants unlike 
from the Sun and SN1987A, and the possibility of dom- 
inant axion emis sion from red giants is not excluded by 
observation [HI Raffelt (2008)]. For the DFSZ axion, 
the region m a > 10~ 2 is excluded due to the large axion- 
electron coupling. For the hadronic axion, Ref. [309l 
Raffelt, Dearborn] argue that the the axion mass greater 
than about 2eV/[(E/N - 1.95)/0.72] would reduce the 
helium-burning time scale and is not allowed. 

For the supernovac explosion, the core temperature 
can go much higher than the temperature in the igni- 
tion phase of Helium in the small mass red giant cores. 
For supernovae, therefore, nuclear reactions are more im- 
portant and the F a bound can be very strong. So we 
use the axion couplings to nucleons discussed in Sub- 
sec. IIIIC II to study the core evolution of supernovae. In 
the beginning, the bounds on the axion decay constant 
were simply obtained by comparing the nucl ear burn- 
ing rates of producing axions and neutrinos J206I . 12891 
Iwamoto (1984), Pantziris, Kang (1986)]. The discov- 
ery of SN1987A was important in that it propelled a 
great deal of i nterest anew in the c alculation of axion pro- 
duction rate [HI SzS EH US Raffelt, Seckel (1988), 
Turner (1988), Mayle, Ellis, Olive, Schramm, Steigman 
(1988), Hatsuda, Yoshimura (1988)]. In principle, the 
same kind of bound on F a could be obtained from the 
earlier supernovae studies. The studies after the dis- 



covery of SN1987A were performed with the derivative 
coupling and quartic terms of Subscc. IIII C II and ob- 
tained a bound F a > 10 9 GeV. But as pointed out in 
IH [MSIS SH, Choi, Kang, Kim (1989), Kang, Pantzi- 
sis (1991), Turner, Kang, Steigman (1989), Carena, Pec- 
cei (1989)] and Subscc. IIII C II a proper treatment of 
nucleon states must be taken into account. For the axion 
emission from supernovae, one must cons train the energy 
output to e„ < 1 x 10 19 crg g"^" 1 [lo3, Raffelt (1990)]. 

Raffelt 



F (86) 



The axion emission rate calculation of Ref. [30 
(2008)] is 



e a = 3.0 x 10 37 [erg g'V 1 ] C 2 N F-*T* 



30MoV 



where F Q , GoV = F a /GcV, T Q , 3 oMcV = T/(30 MeV), and 
F = 0(1). In the supernovae explosion the axion emis- 
sion can be comparable to neutrino emission. Such rem- 
nant axions from all the past supernovae explosions may 
be around us but will be difficult to be detected because 
of the small 1/F a [HI Raffelt (2008)]. For the smaller F a 
region from supernovae explosion, axions can be trapped 
if the axion-nucleon interaction is strong enough. For 
the h adro nic axion, it gives the bound on m a > 1 eV 
[30J HI Turner (1990), Raffelt (1990)], and we have 
the hadronic axion window in the eV range. 

For the KSVZ axion and the MI superstring axion, 
C\ terms are present. For example we can simply take 
5i = 3 and cf = |, corresponding to Z — 0.5, and hence 
obtain c app 



3 F + T1 S 



0.17 for the KSVZ axion. 
Using c app as CV in Eq. (|86j) . we obtain an F a bound 
from supernovae, 



F a > 0.5 x 10 9 GeV. 



(87) 



The white dwarfs in the final evolutionary stage of 
low mass stars (M < 10 ± 2 Mq), with the theoreti- 
cal model implemented in t he DFSZ model, may give a 
stronger bound on F a [303l . Raffelt (1986)] for some re- 
gion of the DFSZ parameter tan/3 = (H u ) / (H^). The re- 
cent study of the bremsstrahlung process gives the bound 
F a > 0.6 x 10 10 sin 2 (3 GeV, and even fits the cooling di- 
agram nicely with F a ~ 1.2 X 10 10 x sin 2 (3 GeV for Hd 
giving mass to electron 205, Iscrn, Garcfa-Bcrro, Torres, 
Catalan (2008)]. Note that tan/3 is known to be large 
(> 30) in SUSY GUT models, and the white dwarfs may 
give the strongest F a bound for some DFSZ axion mod- 
els. 

The axion-nucleon coupling gets enhanced in strong 
magnetic field. Magnetic fields as strong as B > 10 18 
Gauss in neutr on st ars have been assumed in the scalar 
virial theorem [2Zl Woltjer (1964)]. With B > 10 20 
Gauss at the surface, the axion emission rate from neu- 
tron stars or white dwarfs will get enhanced by 0(1) 
compared to the B = case (l97l . Hong (1998)]. 

In summary, axions once produced in the hot plasma 
of a star most probably escape the core, taking out en- 
ergy. This contributes to the energy loss mechanism of 
a star and is used to constrain axion models. From the 
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nucleon-nucleon-a coupling, SN1987A gives the strongest 
astrophysical bound onthe axion decay constant, F a > 
0.5 x 10 9 GeV J2Ql ErE [HI Turner (1990), Raffelt 
(1990), Raffelt (2008)]. 



B. Axions in the universe 



(2005)]. So, with SUSY the domain wall is not so prob- 
lematic. In this case, the problem of strin g radiated ax- 
ions requiring axion mass m a > 10~ 3 eV |ll7l Il2ll . Il9ll 
Davis (1985), Harari, Sikivie (1987), Dabholkar, Quash- 
nock (1990)] is no longer problematic. 



Axions with m a > 24 eV have a lifetime shorter 
than the age of the universe. In this case, axion de- 
cay can lead to photons that can be tested against the 
observed electromagnetic background of the universe, as 
in some spontaneo usly broken flavor symmetric models, 
vi — > Vjd — > ia/77 [39i . Berezhiani, Khlopov, Khomeriki 
(1990)]. However, in this case the needed decay constant, 
10 6 GeV, is outside the current bound on F a . 

The axion for m a < 24 eV has a longer lifetime than 
the age of the universe and can affect its evolution. The 
heavy thermal axions around the eV mass range of Fig. 
[2Kb) become the hot DM in the universe. For the 3- 
8 eV mass range, they accumulate in galaxy clusters 
where their slow decay produces a sharp line that, in 
principle, can be observed by telescope searches as sug- 
gested in [4l|, Bershady, Ressell, Turner (1991)]. In this 
case, the neutrino and axion hot DM must be consid- 
ered together, w hich now constrains the axion mass to 
m a < 1.02 eV [HI, [HI Hannestad, Mirizzi, Raffelt, 
Wong (2007, 2008)], almost closing the hadronic axion 
window of 1-20 eV of Fig. [15j 

But more attention is paid to axions behaving as the 
CDM candidate. The axion potential is almost flat as de- 
picted in Fig. [TT] Therefore, a chosen vacuum stays there 
for a long time, and starts to oscillate when the Hubble 
time H" 1 is comparable to the oscillation period (the 
inverse axion mass), 3H « m a . This occurs when the 
temperature of the universe is about 1 GeV [l|, Il33l 13021 
Preskill, Wise, Wilczek (1983), Abbott, Sikivie (1983), 
Dine, Fischler (1983)]. There exists the doma in wall 
problem in the standard Big Bang cosmology [3321 Sikivie 
(1982)]. The axion strings and the doma in wall problem 
have been elegantly summarized in |336l . Sikivie (2008)]. 
The axion cosmology is correlated to the reheating tem- 
perature Tr,h in the inflationary models, where one must 
deal with both the inflaton and the axion. The density 
perturbations produced by the perturbations of the infla- 
ton field is adiabatic, Sp ma , tter /p ma , ttei = {3/4)6p ia d/p ia d- 
On the other hand, the perturbations produced by the 
fluctuations of the axion field is isocurvature. If the re- 
heating temperature Trh is above the axion scale F a , the 
limit on the isocurvature less than 30% from the large 
scale structure data can be used [38|, Beltran, Garcia- 
Bcllido, Lesgourgues (2007)]. This will be commented 
more in Subsec. IIV CI on the anthropic argument. 

In supersymmetric models, the reheating temperature 
is constrained to T RH < 10 9 GeV or 10 7 GeV (if the 
gluino is lighter than the gravitino) from nucle osynthesis 
requirements in models with a heavy gravitino [1551 . 12171 
Ellis, Kim, Nanopoulos (1984), Kawasaki, Kohri, Moroi 



FIG. 11: The almost flat axion potential. The misalignment 
angle is expected to be of order 1 but can also be very small 
as shown by the thick blue arrow. 

Axions are created at T ~ F a , but the axion vacuum 
(a) does not begin to roll until the Hubble parameter 
reaches the axion mass 3H = m a , which occurs at T ~ 1 
GeV. From then on, the classical field (a) starts to oscil- 
late. For a small misalignment angle, the energy density 
behaves like that in the harmonic oscillator m^F^ which 
is proportional to the axion mass times the number den- 
sity. Thus, its behavior is like that of CDM, which is 
the reason that the axion DM is CDM even though its 
mass is very small and its interaction strength is much 
weaker than "weak" . Even for a large misalignment an- 
gle, an adiabatic invariant / exists and one can estimate 
the current axion energy density. The axion field evo- 
lution with the adiabatic change of the axion mass has 
been considered before [Z3,[73|, Chang, Hagmann, Sikivie 
(1999, 1998)]. 

The temperature dependent axion mass [1831 . Gross, 
Pisarski, Yaffe (1981)] enters in the determination of cos- 
mic temperature T\ where 3H(Ti) ~ m a (Ti). The new 
estimate of T x for F a < 10 16 GeV is a bit below 1 GeV, 
Ti ~ 0.92 GeV [U, Bae, Huh, Kim (2008)]. QCD has two 
phases: the quark-gluon phase and the chiral symmetry 
breaking hadronic phase. Near the critical temperature 
T c , these two phases are separated above and below T c . 
The critical temperature is estimated as 148^^ (172^34) 
MeV for three (two) light quark flavors [52|, Braun, Gies 
(2007)]. So, cosmology near T c needs an information on 
the temperature dependent axion mass. This region is 
in the boundary of weak and strong coupling regimes 
and it is very difficult to estimate the axion mass ac- 
curately. Ear ly attempts in this direction are given in 
[HH HI3, HH, Steinhardt, Turner (1983), Seckel, Turner 
(1985), Turner (1986)]. 

The 't Hooft dctcrminental interaction is given 
schematically in Fig. [SJ In the quark-gluon phase, we 
have the first diagram in the box, which is parametrized 
as —K~ 5 (m u m c im s /pP) cos [( 02+03)6'] where p is the effec- 
tive instanton size in the instanton size integration. Ref. 
PM Gross, Pisarski, Yaffe (1981), Eq. (6.15)] expresses 
the result as 

n( P , ) e {-^ 2 (2N+N f )-12A(\)ll+i(N~N f )}} (88) 

where A = irpT, A(X) ~ -^hi(l + A 2 /3) + a(l + 
7A -2/3-)-8 with a = 0.01289764 and 7 = 0.15858, and 
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the pre-factor n(p, 0) is the zero temperature density 



n(p,0) = m u m d m s C N (£p) 3 — 



4tt 2 V 



Wp)J 



(89) 



Here, the parameters are £ = 1.33876 and CV = 0.097163 
for N = 3. For the QCD coupling constant, we use the 
three loop result @, Amsler et. al. (2008), QCD by 
Hinchliffe], 



a c (p) 



4tt - /3 ln( M VA^ CD ) 



2ft ln[l"(M 2 MqcD)] 



+ 



pfe X ( ( ln [ln(M 2 /A| CD )] - |) 



/3 2 /3o 



(90) 



where & = 11 - |iV/,/3i = 51 - fN f , and /3 2 = 2857 - 
^N f + ^-N 2 f . At T = T GoV GcV (from 700 MeV to 1.3 
GeV), let us parametrize the instanton size integration of 
as 



V(9) = -C(T) cos(0), (91) 

where 6> = a/F a and C(T) is 

C(T) = a inst GeV 4 (T GcV )- n . (92) 
We obtain a inst = 2.858 x 10" 12 (9.184 x 10~ 12 , 7.185 x 
I 



10" 13 ),n = 6.878 
380 (440, 320) MeV [H 



(6.789,6.967) for A QC d = 
Bae, Huh, Kim (2008)]. Equat- 



ing 3H(T) and m{T) = V /C(T)/F 2 we obtain the fol- 
lowing Ti for Aqcd = 380 MeV Hi Bae, Huh, Kim 

(2008)], 



T x 



GcV 



= 0.889 



F„ 



GeV 



10 



12 



(93) 



For F a = 10 12 GeV, we obtain T x ~ 0.89 GcV. This 
number is smaller than those given in 1980s because we 
used the smaller number for the product of current quark 
masses m„m c i rn s b ased on the recent compilation of light 
quark masses [268, Manohar, Sachrajda (2008)]. 

(l)No sudden change of m a (T): 

Since the potential varies much more slowly than the 
field itself, we can use the so-called 'adiabatic invariant 
theorem' that if the potential is adiabatically changed, 
the area in the phase space swept by the peri odic mo- 
tion is unchanged per one axion oscillation [25ll . Landau, 
Lifshitz (1976)]. In this case, for a small misalignment 
angle the adiabatic invariant is p(t)/m(t) which can be 
interpreted as the conservation of the total axion number. 
For a large 9\ , the invariant is not the axion number den- 
sity, but the CDM energy density which can be related 
to the axion number density by a correction factor [19l . 
Bae, Huh, Kim (2008)]. If wc apply this until now, we 
obtain 



Pa (T 7 = 2.73K) = m a (T 7 )n a (T 7 )/i(0 2 ) 



l + Z 



3 ■ 1.66g„(TjT* F a etfti 



2y/g^)M F ^ 



7 



-3-?i/2 



(94) 



r 



where fi(02) is the anharmonic correction and we used 
Z = m u /m d ~ 0.5, m w = 135.5MeV, /„. = 93MeV and 
(^(present) = 3.91. 7 is the entropy increase ratio from 
extra particles beyond the SM. This becomes roughly 

1.449 x 10- 11 £ ( 10 i^'; eV ) F(6 u n) eV 4 , where 9 X is 
the initial misalignment angle at T\ and #2 is the an- 
gle at somewhat lower temperature T2 where the adi- 
abatic invariant / is calculated. The total correction 
factor F{9\,n) takes into account the anharmonic ef- 
fect and the initial overshoot of the misalignment an- 
gle, presented in [H, Bae, Huh, Kim (2008)]. For the 
critical density p c = 3.9784 x 10- n (/i/0.701) 2 (cV) 4 
and Aqcd = 380 =F 60 MeV, the axion energy fraction, 
in terms of F a only, is [l§ Bae, Huh, Kim (2008)] 



n a = 0.397ABC 



nnoi) 



0.701 



(95) 



where A = (m u m d m s /3 ■ 6 • 103 
(F a /10 12 GcV) 1 - 184 - - 010 * with x = 
1, and C = (A QCD /380 MeV) 



-0.733 



MeV 3 )" 0092 , B 
(A Q cd/380 McV) 



(2) Sudden change of m a (T): 

Let us try to calculate the misalignment angle be- 
low the critical temperature of chiral symmetry break- 
ing where a sudden phase change is experienced near the 
critical temperature T c . The QCD interaction for light 
quarks below 1 GeV can be written as 



C = - [m u u L u R + m d d L d R + m s s L s R + h.c. 
~K- b {u L u R d L d R s L s R e- 1 ^ 6 + h.c.) 



(96) 



where K has the mass dimension arising from QCD in- 
stanton physics. The 't Hooft determinental interaction 
[346l 't Hooft (1976)] written above is equivalent to the 
anomaly term and has the same chiral symmetry bchav- 
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ior. For T c < T < F a , quark bilinears are not developing 
VEVs, and the relevant determinental interaction for the 
axion is the first diagram inside the box of Fig. [9l Now, 
the importance of the determination of T\ is how much 
the misalignment angle 9\ can be shrunk at T c . 



V 
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m l« 
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FIG. 12: Phase transition near the critical temperature T c w 
150 MeV. 

In the hadronic phase below the critical temperature 
T c . the axion potential is schematically depicted in Fig. 
IT2l The value 0i(T c ) is the boundary value at T c we use 
in the effective Lagrangian below T c . Below T c , the quark 
bilinears develop VEVs and we must consider the possi- 
bilities of qLlFt. replaced with ((Jl<1fi)- The effective La- 
grangian from the determinental interaction is depicted 
in Fig. M 

In the limit F a 3> /', the mass eigenstatcs in one flavor 
QCD arc 



1 

f'/Fg 

l+m/K' 



-f'/Fa 
l+m/K' 

1 



(97) 



Eq. (|38|) with v 3 — has minima at 9 = 2irn (n = 
integer). For v 3 ^ 0, minima are at 9 V ' = 2-Km (m = 
integer) and 9 = 2ttti (n = integer). Therefore, the 9 
direction can be taken as the approximate axion direction 
even below T c . The minimum point in the direction of the 
axion is not changed when one goes from 9^ ^ to 9 V > = 
0, i.e. above and below the critical temperature. [If the 
minimum of 9 were shifted by ir from going from 9 n i = 
to 6 v i = 7r, the shrunk 9\{T C ) at T c is near ir, and we 
must start from 0(1) misalignment angle at T c .\ In most 
regions of the phase transition space, a time scale At is 
needed for the sound wave of quark bilinears to propagate 
to a large distance, which releases the latent heat to keep 
the constant tem perature during the first order phase 
transition in |277l . Mukhanov (2005), p. 149]. Even if one 
considers supercooling toward a sudden phase transition, 
the parameter space for a sudden phase change is almost 
nill and the axion energy density presented in ([93)1 is 
reliable [H, Bae, Huh, Kim (2008)]. 

In Fig. 1131 we present the exclusion plot for m u = 
2.55 MeV, m d = 5.04 MeV and m s = 104 MeV [26J 
Manohar, Sachrajda (2008)] in the F a versus 9\j ^H 



FIG. 13: The plot for F a versus the misalignment angle Oi/y/y 
as a function of il a . The overclosure p ortio n is from the preci- 
sion measurement requiring Q a < 0.23 [2471 ]. The green region 
is the excluded region by the condition Q a > 0.23. The yellow 
band is the error bar region of Aqcd and two red lines are 
the limits from the light quark mass bounds. 



space, incl uding the anharmonic effect and the WMAP 
value [HI Dunkley et al. (WMAP Collaboration, 2008)] 
of CDM density combined with additional data [247l . Ko- 
matsu et al. (WMAP Collaboration, 2008)] f! DM /i 2 ~ 
0.1143 ± 0.0034. Note that F a of order 10 13 GeV is 
not very unnatural, whi ch results from the new smaller 
masses for u and d [268l Manohar, Sachrajda (2008)]. 



If the axion is the CDM component of the universe, 
then they can be detected even though it may be very 
difficult. The feeble axion coupling can be compensated 
by the huge number of axions, since the number density 
is ~ F% and the cross section is ~ 1/-Ff. So, there is 
a hope to detect cosmic axio ns, w hich has been realized 
by Sikivie's cavity detector [333l Sikivie (1983)]. But 
the Sikivie detector has the technical limitation for the 
interesting large and median regions of the F a window. 
For example, the F a region F a > 10 13 GeV advocated 
in anthropic arguments needs too large cavity size and 
the supergravity mediation preferred region F a ~ 5 x 
10 10 GeV requires 1.6 mm order cavities. For technically 
preferred axion masses in the region 10 -6 eV, one needs 
a low temperature cavity with dimension 0{> 10 4 cm 3 ) 
and a magnetic field strength of O(10 Tcsla). The current 
status of cosmic axion search is shown in Fig. 1141 
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FIG. 14: The bounds on cosmic axion searches with some 
theoretical predictions. The coupling on the vertical axis is 
the coefficient of E • B. The future CARRACK and ADMX 
experiments are from J20ll . |354| . 



C. Axion cosmology beyond the window 

If F a > 10 12 GcV, an 0(1) misalignment angle Q\ 
is ruled out from the cosmic energy density argument. 
However, if ^i -C 1, the axion energy density can be 
within the closure density. Rather th an fi ne-tuning 6\ 
to order 10~ 3 for a Planck scale F a [293, Pi (1984)], 



the anthropic argument of Weinberg [263, 363, Weinberg 
(1987), Linde (1988)], that life forms can evolve in a uni- 
verse with a sufficiently long lifetime, can be used for an 
allowable 6\. 

The homogeneous axion field value (with a — ► —a sym- 
metry) right after inflation can take any value between 
and nF a or 6\ = [0, ir] because the height of the axion po- 



tential is negligible compared to the total energy density 
right after inflation. So, in the axion context with only 
the misalignment production of axions, the CDM den- 
sity is chosen as a random number by the spontaneous 
symmetry breaking of the U(1)pq. Even in the multi 
component CDM models, including axion, the axion mis- 
alignment angle can play as the random num ber . This 
singles out axion physics, as stressed in [345, Tcgmark, 
Aguirre, Rees, Wilczek (2006)], from all other anthropic 
arguments without axion in that selecting an axion vac- 
uum has an unavoidable random process that fixes the 
key cosmological parameter. This also distinguishes ax- 
ions from WIMPs, super- WIMPs etc. where the abun- 
dance is fixed by particle physics parameters and not by a 
primordial random process. So £l a may be at the required 
value by an appropriate initial misalignment angle in 
models with axions with F a > 10 12 GeV. Tcgmark et al. 
studied the landscape scenario for 31 dimcnsionless pa- 
rameters and some dimensionful parameters with which 
habitable planets are considered for the assumed nuclear 
physics paramet ers kiCl Barr, Seckel (1992)]. For exam- 
ple, Fig. 12 of [3451 . Tegmark, Aguirre, Rees, Wilczek 
(2006)] presents the scalar fluctuation Q ~ Sp/p vs the 
matter density per CMB photon £, in which the anthrop- 
ically chosen point is shown as the star. In models with 
axion, this point results from a random number after in- 
flation. If a WIMP is the sole candidate for CDM, one 
obtains just one number for Sp/p from particle physics 
parameters, which may not fit to the observed point of 
that figure. So, we may need the CDM favored WIMP 
and in addition the axion with F a > 10 12 GeV, with the 
axion CDM fraction R a = f2 a /OcDM- But this large F a 
anthropic region has a potential conflict with the WMAP 
5 year data, as presented in the F a v s. Ej(= the inflation 
energy scale) plane of Fig. 2 of [l94l . Herzburg, Tegmark, 
Wilczek (2008)]. For R a = 1, for example, F a > 10 14 
GeV is inconsistent with the WMAP 5 year data on the 
uppe r bound on the isocurvature fluctuation a a < 0.072 
J24J Komatsu et al. (2008)]. 

From the study of outer space axions, we present a 
cartoon for the F a bound in Fig. [T5] where the future 
CAST and ADMX experimental regions are also marked. 



r 



D. Quintessential axion 



In light of SUSY breaking in supergravity, it is gen- 
erally believed that at least a hidden confining force is 
needed at an intermediate scale. This hidden sector and 
the observable sector couple extremely weakly in most 
phenomcnological models. This scheme fits very well in 
the heterotic string framework and in heterotic M-theory. 
In cosmology, on the other hand, we have the i mpo rtant 
dark energy problem already for a decade [2961 . 13151 . Per- 
mutter et al. (1998), Riess et al. (1998)], which has led 



to a lo t of interest in quintessence models since the late 
1980s HH Wetterich (1988)]. The quintessence related 
to axion physics is cal led "quintessential axion" (Q A) 
which was suggested in [24 ll . Kim, Nilles (2003)]. There 
were attempts to identi fy on e of the MD axions as the 
quintessential axion [13, [IH Kim (2000), Choi (2000)]. 

To explain the dark energy in terms of a QA, one re- 
quires the VEV of QA not to roll down until recently. Of 
course, it is required for the current vacuum energy den- 
sity of the classical QA to be of order A 4 w (0.003 eV) 4 . 
These two conditions restrict the QA decay constant f q 
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FIG. 15: A cartoon for the F a bounds. 



and the QA mass m q . 
potential as 



We can parametrize the QA (<j>) 



For 



v[4>] = A 4 t/(o, e 



P/P 



fa 



(98) 



< —1 + 5, we require f q > 
v/(2 - 6)/ 66 M P \U'\ where U' = dU/d£ [HE Kim, 
Miles (2003)]. Generically, one needs a Planckian scale 
quintessential axion decay constant f q . So, the QA mass 
is extremely small, < IO -32 eV. As a result, there are two 
problems to be resolved to achieve the QA idea: a large 
decay constant and an extremely shallow QA potential. where c, 



V. AXION DETECTION EXPERIMENTS 

There are currently a variety of experiments searching 
for axions, whether they are left over from the big bang 
or produced in stars or the laboratory. Though these 
experiments search for axions at a variety of mass and 
coupling scales they all rely on the Primakoff process for 
which the following coupling, c a77 is given in Eq. (|75[1 . 



C 



Fa 



1.98 



(101) 
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Solar axion search 



It has long been believed that the MI axion has rather 
a robust mod el independent prediction of its decay con- 
stant HI,!!!! Choi, Kim (1985), Svrcek, Wittcn (2006)]. 
Recently, however, it was shown that the MI axion may 
not be model independent since the decay constant may 
depend on the compactification scheme in warped inter- 
nal space, ds 2 = h 2 u r] fjtu dx ,J, dx 1 ' + g m n(y)dy m dy n [ll9L 
Dasgupta, Firouzjahi, Gwyn (2008)], 



f - H™L 

a " V M P 



(99) 



where /3 depends on the warping in the compact space 
y£K, 







J dPy^/g^hl 



(100) 



Thus, the MI axion with a small (3 can be a QA if the 
QCD axion decay constant can be in the intermediate 
scale. This possibility may be realizable in some com- 
posite axion models as recently suggested in [2421 . Kim, 
Nilles (2009)]. 



1. Axion Helioscopes 

Axions produced in the nuclear core of the sun will 
free-stream out and can possibly be detected on Earth 
via an axion helioscope, first described in 1983 |333l . l334L 
Sikivie (1983, 1985)] an d dev eloped into a practical labo- 
ratory detector in 1988 [355l . van Bibber, Mclntyre, Mor- 
ris, Raffelt (1989)]. The technique relies on solar axions 
converting into low energy X-rays as they pass through 
a strong magnetic field. The flux of axions produced in 
the sun is expected to follow a thermal distribution with 
a mean energy of (E) = 4.2 keV. The integrated flux at 
Earth is expected to be <& a = q 2 n 3.6 7 x 1 11 cm _2 s _1 
with <7i = (a C m/27rF a )c a77 10 10 GeV J38CJ Ziotas et al. 
(2005)] . The probability of a solar axion converting into a 
photon as it passes through a magnet with field strength 
B and length L is given as: 



P = 



^cm^a77 -^ ^ 
4^Fa 



2L- 



1 — cos(qL) 

{qL) 2 ' 



(102) 



Here c a77 is defined as the coupling of the axion to two 
photons as given in Eq. (|101[) . while q is the momentum 
difference between the axion and the photon, defined as 
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q = m\/2E where E is the photon energy. To maintain 
maximum conversion probability the axion and photon 
fields need to remain in phas e ove r the length of the mag- 
net, thus requiring qL < tt [355, van Bibber, Mclntyre, 
Morris, Raffelt (1989)]. For low mass axions q — > lead- 
ing to a maximum conversion probability. More massive 
axions will begin to move out of phase with the photon 
waves though this can be compensated for by the additon 
of a buffer gas to the magnet volume, t hus i mparting an 
effective mass to the conversion photon |355l van Bibber, 
Mclntyre, Morris, Raffelt (1989)] and bringing the con- 
version probability back to the maximum. Various axion 
masses can be tuned to by varying the gas pressure. 

An initial axion helioscope was built at Brookhaven in 
1992 and used a 2.2 ton iron core dipole magnet oriented 
at Su n with a proportional chamber for X-ray detection 
[2561 . Lazarus et al. (RBF Collaboration, 1992)]. It was 
followed by a 4 tesla superconducting helioscope, devel- 
oped by the University of Tok yo, which ran for 1 week 
with an evacuated bore in 1997 275l . l287l . Moriyama et al. 
(1998), Ootani et al. (1999)] and for 1 month with a he- 
lium filled bore in 2000 ^M Inoue et al. (2002)]. Though 
both managed to set limits over a wide mass range their 
sensitivities were still well above even the most optimistic 
KSVZ axion couplings. Recently though, the University 
of Tokyo group released data taken between December 
2007 and April 2008, which was able to set a limit of 



3 a77 < 5.6 - 13.4 x 10 
mass range 0. 



10 GcV 1 f or th e axion in the 
C m a < 1.00 eV [lol Inoue et. al. 
(Tokyo Axion Helioscope Collaboration, 2008)]. 

In order to push into proposed axion model space 
third generation axion helioscopes have been developed 
at CERN (the CERN Axion Solar Telescope (CAST)) 
and at the University of Tokyo. Utilizing a prototype 
Large Hadron Collider (LHC) magnet with L = 9.3 m 
and B = 9 tesla CAST began taking data in 2003. It 
utilizes a rail system to track Sun for 90 minutes a day 
at sunrise and sunset and its dual magnet bore allows 
it to employ up to four different X-ray detectors (one 
on each end of each magnet bore). Currently a time- 
projection chamber (TPC), a Micromegas (micromesh 
gaseous structure) detector and an X-ray reflective tele- 
scope with a charge coupled device (CCD) detector are all 
used to detect converted X-rays. Results from the com- 
bined 2003 and 2004 runs yield limits on axion-photon- 
photon couplings to c ail /F a < 7.6 x 10~ 8 GcV _1 (lol 
Andriamonje et al. (CAST Collaboration, 2007)]. The 
experiments second phase utilizing 4 He and 3 Hc buffer 
gases is currently underway with the latter gas allowing 
for axion searches in proposed model space up to mass ~ 
1 eV. 



the Bragg conditons [29Gt Paschos, Zioutas (1994)]. Var- 
ious dark matter WIMP search collaborations were able 
to look through their data sets and set limits on possible 
interactions from solar axions. Th ese in cluded germa- 
nium experiments such as COSME |274l . Morales et al. 
(2002)] and SOLAX 0, Avignone et al. (1998)], the 
reactor germanium experiment TEXONO ITOL Chang et 
al.], as well as the DAMA experiment [40l.l42l. Bernabei 
et al. (2001, 2003)] which utilized Nal crystals. The lim- 
its from these searches can all be seen in Fig. [TBI One 
advantage of this technique is that their sensitivity is in- 
dependent of axion mass, as long as one can neglect any 
nuclear recoils [6fJ, Carosi, van Bibber (2007)]. 

Current results for solar axion searches can be seen in 
Fig. M 
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FIG. 16: Exclusion plot of axion-photon coupling versus axion 
mass [6lj . The black bold line limit is for the phase one of the 
CAST experiment and results with inclusion of buffer gas are 
expected to push up in mass reach to plausible axion models. 
The field theoretic expectations are shown together with the 
string theory Zrz-i model of [94| . In the lower left apricot 
box, Fig. [14] is located. 



2. Bragg Diffraction Scattering 



3. Geomagnetic Conversion 



An alternative to axion helioscopes was proposed in 
1994, using crystal detectors to search for X-rays gener- 
ated by coherent axion-to-photon conversion which meet 



It has recently been pointed out J122L Davoudiasl, Hu- 
ber (2006)] that solar axions might pass through Earth 
and convert to X-rays on the other side as they pass 
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FIG. 17: Outline of the general configuration of a resonant 
microwave cavity detector along with the associated singal ex- 
pected from axion-photon conversions. This diagram includes 
both the virial component as well as possible lines from co- 
herent axions. 



through the earth's magnetic field. They could then be 
detected by X-ray telescopes and the solar X-ray back- 
ground could be effectively shielded by Earth. 



B. Search for cosmic axions 

Cosmic axions left over from the Big Bang may be de- 
tected utilizing microwave cavity haloscopes [3331 . I334L 
Sikivie (1983, 1985)]. The strategy relies on primordial 
axions drifting through a microwave cavity immersed in a 
strong static magnetic field in which they can resonantly 
convert to microwave photons. The cosmic axions feeble 
interactions can be in part compensated by their large 
numbers; since the number density goes as ~ F^ while 
their cross section goes as ~ 1/f a • ^ the axion makes up 
the majority of CDM in the universe, its l ocal density is 
expected to be roughly 0.45 GeV/cm 3 163, Gates, Gyuk, 
Turner (1995)], which yields a number density of ~ 10 14 
axions/cm 3 if one assumes a 4.5 pcY axion. The expected 
microwave signal will be a quasi-monochromatic line be- 
ginning at the microwave frequency corresponding to the 
axion mass and slightly broadened upward due to the ax- 
ion virial distribution, with expected velocities of order 
10 _3 c implying a spread in energies of SE/E ~ 10~ 6 . 

There could also be an additional signal from non- 
thermalized axions falling into the galaxy's gravitational 
well which would yield very sharp signals due to their 
low predicted velocity dispersion (< 10 _7 c) [335l Sikivie 
(2003)]. 



1. General Detector Properties 

Since the Lagrangian for axions coupling to a magnetic 
field goes as 
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the only resonant modes which can couple to axions 
are those that provide an axial electric field component 
(transverse magnetic or TM modes). The expected power 
generated f rom axion-to-photon conversions in the cavity 
is given by [334l . Sikivie (1985)] 
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where V is the cavity volume, Bq is the magnetic field 
strength, p is the local axion mass density, m a is the axion 
mass, Cimn is a form factor which describes the overlap 
of the axial electric and magnetic fields of a particular 
TMi mn mode, Ql is the microwave cavities loaded qual- 
ity factor (defined as center frequency over bandwidth) 
and Q a is axion quality factor defined as the axion mass 
over the axion's kinetic energy spread. The mode depen- 
dent cavity form factor is defined as 
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where E^[x)e %UJt is the oscillating electric field of the 
TM\ mn mode, B (x) is the static magnetic field and e 
is the dielectric constant of the cavity space. For a cylin- 
drical cavity with a homogeneous longitudinal B-ficld the 
To i o mode yields the largest form factor with Coio ~ 0.69 
g3, Bradley (2003)]. 

The mass range of cosmological axions is currently 
constrained between peV and meV scales which corre- 
sponds to converted photon frequencies between several 
hundred MHz and several hundred GHz. Since larger mi- 
crowave cavities correspond to lower resonant frequencies 
and lighter axions are more likely to contribute to the 
dark matter density experiments have been designed to 
start searching at the low end of the frequency range. At 
these frequencies cavities can only scan a few kHz at a 
time in order to maintain the maximum quality factor. 
Axial metallic and dielectric tuning rods are utilized to 
tune the cavities resonant frequency as it scans over the 
possible axion mass range. The scan rate is determined 
by the amount of time it takes for a possible axion signal 
to be detected over the microwave cavity's intrinsic noise 
and is governed by the Dicke radiometer equation [1281 . 
Dicke (1946)], 
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Here P a is the power generated by axion-photon conver- 
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sions (Eq. (|104|) ). Pn = ksBTs is the cavity noise power, 
B is the signal bandwidth, t is the integration time, kg is 
Boltzmann's constant and T$ is the system temperature 
(electronic plus physical temperature) . The scan rate for 
a given signal to noise is given by 
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One can see from Eq. (|106p that even a tiny expected 
signal power can be made detectable by either increasing 
the signal power (P a ex VBq), increasing the integration 
time t or minimizing the system noise temperature T$- 
Technology and costs limit the size and strength of the 
external magnets and cavities and integration times are 
usually t ~ 100 seconds in order to scan an appreciable 
bandwidth in a reasonable amount of time. As a result 
the majority of development has focused on lowering the 
intrinsic noise of the first stage cyrogenic amplifiers. 
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FIG. 18: Schematic of the ADMX experiment [6C 



2. Microwave Receiver Detectors 

Initial experiments wer e undertaken at Brookhaven 
National Laboratory |l26t DePanfilis et al. (Rochester- 
Brookhaven Collaboration, 1987)] and the University of 
Florida [HI Hagmann et al. (U. of Florida, 1990)], but 
their modest sized cavities and magnet fields meant they 
were still factors of 10-100 times away from plausible ax- 
ion model space. There are currenly two active second 
generation experiments underway, the Axion Dark Mat- 
ter eXpereiment (ADMX) at Lawrence Livermore Na- 
tional Laboratory (LLNL), USA and the Cosmic Ax- 
ion Research with Rydberg Atoms in Cavities at Kyoto 
(CARRACK) experiment in Japan. Both experiments 
utilize large microwave cavities immersed in a strong 
static magnetic field to resonantly convert axions to pho- 
tons but they go about detecting these photons in two 
different ways. ADMX uses ultra-sensitive microwave re- 
ceivers while CARRACK uses Rydberg atoms to detect 
single photons. 

The ADMX experiment is a collaboration of LLNL, 
MIT, the University of Florida, Lawrence Berkeley Na- 
tional Laboratory (LBNL), U.C. Berkeley, U. of Chicago 
and Fcrmilab, and has been operating in various modes 
since February, 1996. A diagram of the experiment can 
be seen in Fig. [T5] ADMX consists of an 8.5 Tesla super- 
conducting magnet, 110 cm in length with a 60 cm clear 
bore. A 200 liter stainless steel microwave cavity plated 
in ultra-pure copper is suspended below a cryogenic stage 
in the center of the B-field. Power generated in the cav- 
ity is coupled to an adjustable antenna vertically input 



through the top cavity plate. Any signal is then boosted 
by extremely low noise cryogenic amplifiers before being 
sent through a double-heterodyne mixing stage. Here the 
GHz range signal is mixed down to an intermediate 10.7 
MHz, sent through a crystal bandpass filter, and then 
mixed down to audio frequencies at 35 kHz. This audio 
signal is then analyzed by fast-Fourier-transform (FFT) 
electronics which measure over a 50 kHz bandwidth cen- 
tered at 35 kHz. There is also a "high resolution" chan- 
nel in which the signal is mixed down to 5 kHz and sent 
through a 6 kHz wide bandpass filter. Time traces of the 
voltage output, consisting of 2 20 data points, taken with 
a sampling frequency of 20 kHz is then taken, resu lting 
in a 52.4 second sample with 0.019 Hz resolution J149L 
Duffy et al. (2006)]. 

Since the system noise is dominated by the first stage of 
amplification great care was taken in choosing the cryo- 
genic amplifiers. The initial ADMX data runs utilized 
Hetero junction Field Effect Transistor (HFET) ampli- 
fiers developed by th e National Radio Astronomy Obser- 
vatory (NRAO) [lH, Daw, Bradley (1997)]. Even though 
they had noise temperatures of only 2 K, the quantum 
noise limit at a GHz (defined as T q = hv/ks) is only 
50 mK. As a result a great deal of development went into 
replacing the HFETs with more sensitive Superconduct- 
ing Quantum Interference Devices (SQuIDs) which had 
noise temperatures of only 15% the quantum limit [471 . 
Bradley (2003)]. As of this writing data is being taken 
using the SQuIDs for the first stage of amplification. 

Results from the initial run using HFET amplifiers 
have already probed plausible axion model space in the 
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FIG. 19: Medium resolution limits from the ADMX experi- 
ment. The left-side plot shows limits to axion coupling assum- 
ing a dark matter halo density of p = 0.45 GeV/cm 3 (upper 
region excluded). The right-side plot shows limits on the ax- 
ion contribution to the dark matter halo density assuming the 
axion has either the KSVZ or the DFSZ couplings [47J . 
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FIG. 20: High resolution limits from the ADMX experiment. 
Limits given in terms of axion contribution to the dark matter 
halo density assuming the axion has either the KSVZ or the 
DFSZ coupling strength. The medium resolution plot fo r tha t 
mass range for the DFSZ coupling is also given as red [l49( |. 
The KSVZ axion with eg = shown above gives p a at Earth 
less than 0.16 GeV/cm 3 which corresponds to Q a < 0.36. So 
the ADMX line of Fig. [14] using Eq. ([95]) crosses the e Q - 
KSVZ line and goes down to the DFSZ line. 



axion mass range between 2.3 - 3.4 fi eV [471 . Bradley 
(2003)]. Results from a high resolution search have 
probed further into c oupli ng space over a smaller mass 
range, 1.98-2.17 AieV [ill Duffy et al. (2006)]. As of 
this review ADMX is scanning over the mass range cor- 
responding to 800-900 MHz using SQuID amplifiers. 



3. Rydberg Atom Detectors 

The CARRACK experiment has published proof of 
concept papers for their detection technique using Ryd- 
berg atoms as opposed to low noise amplifiers [344], Tada 
et al. (2006)]. The experimental setup can be seen in 
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FIG. 21: General schematic of CARRACK experiment utiliz- 
ing Rydberg atoms to reco ver single photons generated in the 
microwave cavity |60l . 13441 ] . 



Fig. [21] In it rubidium atoms are excited into a Ryd- 
berg state (|0 >— > | n >), and move through a detection 
cavity coupled to an axion conversion cavity. The spac- 
ing between energy levels is tuned to the appropriate fre- 
quency utilizing the Stark effect and the Rydberg atoms 
large dipole transition moment ensures efficient photon 
detection (one photon per atom, \ns >—* \np >). The 
atoms are then subjected to a selective field ionization 
allowing the atoms with the higher energy state (\np >) 
to be detected [6fj, Carosi, van Bibber (2007)]. The ad- 
vantage of this system is that Rydberg atoms act as single 
photon detectors and thus do not suffer from quantum 
noise limitations. 

Though still in the development phase CARRACK has 
already gone through two iterations, CARRACK 1 and 
CARRACK 2, and it has measured cavity emission at 
2527 MHz down to a temperature of 67 mK, which is 
a factor of two below the quantum noise floor for that 
frequency. Their eventual goal is 10 mK [344L Tada et 
al. (2006)]. One disadvantage of this technique is that 
one cannot detect signals finer than the bandpass of the 
cavity, of order ~ 10~ 5 , which negates searches for late 
infall coherent axions. 



C. Laser searches 

In addition to cosmological and solar axion searches 
there is also a class of laboratory axion searches which 
utilize laser photons (7i aS cr) traversing a magnetic field. 
Here the polarized laser photons can scatter off virtual 
photons (7v) provided by the magnetic field and con- 
vert into axions 7i aS er + "f v —>■ a. Currently laser axion 
searches have fallen into two general categories. The first 
technique looks for magneto-optical effects of the vac- 
uum due to polarized laser photons disappearing from 
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the beam as they are converted into axions. The second 
looks for photons converting into axions in the presence 
of a magnetic field, which are then transmitted through 
a wall and converted back into photons by a magnetic 
field on the other side, so called "light shining through 
walls" experiments. 
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There can be the axion-photon-photon anomalous cou- 
pling of the form aE-B. The laser induced axion- like par- 
ticle search employing this coupling has been performed 
since early 1990s by the Rochester-Brookhaven-Fermilab- 
Trieste (RBFT) group [H, Cameron et al. (1991)]. A few 
years ago, the same type experiment by the PVLAS col- 
laboration was perf orme d with an initial positive signal 
with F a ~ 10 6 GeV [Izl Zavattini et al. (2006)] we dis- 
cussed earlier. This has led to some exotic models where 
a vacuum dichroism is achieved by producing axion-likc 
particles as shown in Fig. 1227 a). Because of the non- 
renormalizable interaction implied in Fig. I22f a). one may 
reconcile this model with the astrophysical bound [272l 
Mohapatra, Nasri (2007)]. Or if light milli-charged par- 
ticles are produced in the strong magnetic field as shown 
in Fig. 122(b) a vacuum dichroism is achieved as dis- 
cussed in [111 [HI, [Hi, Gies, Jaeckel, Ringwald (2006), 
Masso, Redondo (2006), Kim (2007)]. Here, the polar- 
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FIG. 22: Possible processes leading to a vacuum dichroism. 

ization of the laser beam is looked for. With more data 
accumulation, there is no convincing evidence for an ax- 
ionlikc particle with F a ~ 10 6 GeV at present , contrary 
to an earlier confusion [7y, llOOl Il5ll . 13751 . 1 3781 . Zavattini 
et al. (2007), Chou et al. (2008), Yoo (2007), Chen et 
al. (2006), Dupays et al. (2005)]. But, this incident has 
lead to the current search of axionlike particles at DESY 
[3161 . Ringwald (for the 4th Petras Workshop) (2008)]. 



2. Light shining through walls 

The "Light shinging through walls" techniqu e for 
searching for axions was first proposed in 1987 by [356l 
van Bibber et al. (1987)] and recently a model study 
has been presented [j, Adler, Gamboa, Mendez, Lopez- 
Sarrion (2008)]. The general experimental layout can 
be seen in Fig. [23] where polarized laser photons pass 
through the magnetic field with E\\B and any converted 
axions (or other psucdoscalar particles) can continue 



FIG. 23: Schemtic of "light shining through walls" experiment 

a. 



through an absorber to be re-converted to photons on 
the other side. 

The probability for a photon to convert into an axion 
as it traverses the "Axion Source" region is given by 
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This is the same probability for an axion to convert back 
into a detectable photon in the "Axion Detector" region 
on the otherside of an absorbor, which leaves the total 
probability for detecting a photon-axion-photon conver- 
sion as P 7 ^ a ^ 7 = Pif (ignoring photon detection ef- 
ficiencies of course) [32, Battesti et al. (2007)]. There 
is a maximum detectable axion mass for these laser ex- 
periments as a result of the oscillation length becoming 
shorter than the magnetic field length causing a degra- 
dation of the form factor F(q) = 1 — cos(qL) / (qL) 2 but 
this can be compensated for by multiple discrete dipoles. 

The first experiment using this technique was per- 
formed by a collaboration consisting of Rochester- 
Brookhaven-Fermilab- Trieste (RBFT) in early 1990s (57l 
Cameron et. al. (1993)]. Using two superconducting 
dipolc magnets (L = 4.4 m and B = 3.7 T) and a laser 
(A = 514 nm and P = 3 W) with an optical cavity pro- 
viding ~ 200 reflections in the axion generating region 
they were able to set upper limits on axion couplings of 
g ail < 6.7 x 10~ 7 GeV (95% C.L.) for pseudoscalars with 
a maximum mass of m a < 10 -3 eV [571 . Cameron et al. 
(1993)]. 

Recent photon regeneration exper iments include the 
BMV collaboration at LULI [317, Robilliard et al. 
(2007)] which uses a short pulsed-fiel d ma gnet and the 
GammcV collaboration at Fcrmilab [lOOl . Chou et al. 
(2008)] which uses a Tevatron dipole magnet (L = 6 m 
and B = 5 T) with an optical barrier in the middle. Both 
of which have ruled out the signal reported by the PVLAS 
(see next section). Fig. |2~4"1 illustrates the current bounds 
from these latest regeneration experiments. Recently it 
has been shown that photon regeneration experiments 
can be resonantly enhanced by encompassing both the 
production and reconver sion magnets in matched Fabry- 
Perot optical resonators [3371 Sikivie, Tanner, van Bibber 
(2007)]. 
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FIG. 24: Curr ent limits on axion coupling from the GammeV 
collaboration [lOOl . |375], 
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FIG. 25: The upper figure shows the effect of dichroism as 
photons converting into axions cause a rotation of the linear 
beams polarization vector by an amount e. The lower figure 
shows virtual axions inducing birefringence in which the linear 
beam acquires ellipticity ty |32l |. 



3. Magneto-Optical Vacuum Effects 

An alternative to the "shining light through walls" 
technique is to look for the indirect effect of photons in 
polarized laser light converting into axions as the beam 
traverses a magnetic field. Fig. [25] illustrates the two 
different ways in which axion interactions can modify 
a polarized laser beam, induced dichroism and vacuum 
birefringence. Vacuum dichroism occurs when a polar- 
ized laser beam passes through a dipolc magnet with the 



electric field component £ at a non-zero angle <j> rela- 
tive to the B. The photon component parallel to the 
B will have a small probability to convert into axions 
causing the polarization vector to rotate by an angle e. 
Vacuum birefringence is due to the induced ellipticity of 
the beam (^) as a result of virtual axions. ft should be 
noted that higher order QED diagrams, or "light-by-light 
scattering" diagrams, are expected to contribute to vac- 
uum birefringence as well. Each of these effects can be 
estimated as 
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in the limit that m^L/Au; -C 1. Here L is the effective 
path lengths, N is the number of paths the light travels in 
the magnetic field, m a is the axion mass, uj is the photon 
energy and 6 is the photon polarization relative to the 
magnetic field [H, Battesti et al. (2007)]. 

The initial experiment looking for magneto-optical vac- 
uum effects was ca rried out by the RBFT Collaboration 
in the early 1990s [H§, Semertzidis et. al. (1990)]. This 
experiment used a single-pass length 8.8 m magnet with 
a magnetic field of B ~ 2.1 T and N = 500. It set a limit 
on the polarization rotation of e < 3.5 x 10~ 10 which was 
still 3 orders of magnitude higher than that expected by 
light-by-light scattering and almost 15 orders of magni- 
tude greater than a m a ~ 10~ 3 eV axion. 

Recently the early PVLAS collaboration reported the 
positive detection of vacuum dichroism. This experiment 
consists of a 1 m long 5 T superconducting magnet with 
a angular frequency O mag of the magnet rotation and a 
6.4 m long Fabry-Perot cavity giving the pass number 
N = 2f2 mag /7r ~ 44, 000. It registered a polarization 
shift of 



e = (3.9 ± 0.5) x 10 _12 rad pass 



(111) 



which translates to an allowed mass range of a neu- 
tral pseudoscalar boson of 1 meV< mj <1.5 meV and 
a coupling streng th o f 1.5 x 10~ 3 GcV -1 < c ail /F a < 
8.6 x 10- 3 GcV _1 Hzl, Zavattini et al. (PVLAS Collabo- 
ration) (2006)]. Tho ugh the r eport of this positive signal 
has been retracted [3771 I378J Zavattini et al. (PVLAS 
Collaboration) (2007, 2008)], the interest it raised has 
lead to a number of more advanced experimental searches 
such as some of the new laser regeneration experiments 
mentioned previously. 
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VI. THEORIES FOR VERY LIGHT AXIONS 



Axions from 



Order of F a 



String theory 

M-theory 

Large extra (n) dimension 
Composite models 

Rcnormalizablc theories 



String scale or Planck scale 

String or the scale of the 11th dimension 

Combination of the fundamental mass Md and extra dimension radius R 

Compositeness scale 

The U(1)pq global symmetry breaking scale 



TABLE I: Natural scales of F a . For n large extra dimension, the Planck mass is Mp ~ Md(R/Md)" ■ 

I 



The axion couplings come in three types: the PQ sym- 
metry preserving derivative coupling c\ term, the PQ 
symmetric C2 term, and the anomalous C3 term. The 
PQ symmetry gives a gluon anomaly and ci + C3 must 
be nonzero. Generally, we can therefore define a as a 
pscudoscalar field without potential terms except the one 
arising from the gluon anomaly under certain basis (for 
example in the C2 = basis), 
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Then, we note that this kind of nonrenormalizablc 
anomalous term can arise in several ways. The natu- 
ral scales of F a are shown in Table |TJ In any case, the 
essence of the axion solution (wherever it originates in 
Table HJ is that the axion VEV (a) seeks 9 — what- 
ever happened before. In this sense it is a cosmological 
solution. The potential arising from the anomaly term 
after integrating out the gluon field is the axion potential 
(with c-2 + c 3 ) sketched in Fig. [7] 



A. SM singlets without SUSY 

A complex SM singlet carrying the PQ charge can ap- 
pear in man y ext ensions of the SM: in grand unified theo- 
ries(GUTs) J3p Wise, Georgi, Glashow (1981)], in com- 
posite models [Hill, [IM Kim (1985), Choi, Kim (1985), 
Babu, Choi, Pa ti, Zh ang (1994)], and in models with ex- 
tra dimensions |l3li Di Leila, Pilaftsis, Raffelt, Zioutas 
(2000)]. 

In the SU(5) G UT, the axion can be embedded in a 
complex 24 = E [HI Wise, Georgi, Glashow (1981)], 
in which case the VEV (S) breaking SU(5) down to the 
SM and hence the axion decay constant is the GUT scale 
and is outside the axion window. On the other hand, 
a complex GUT singlet, whose VEV is not related to 
the GUT scale, can house the axion within the axion 



window. The SUSY generalization of the SU(5) GUT 
axion has been shown to be possible [284l . Nilles, Raby 
(1982)]. Recently, the electrophilic axion in view of the 
white dwarf evolution [2051 . Isern, Garcfa-Bcrro, To rres , 
Catalan (2008)] with the two dark matter scenario J200l . 
Huh, Kim, Kyae (2008)] has been suggested in a SUSY 
flipped SU(5) [H Bae, Huh, Kim, Kyae, Viollier (2008)]. 



B. Composite axions 

A SM singlet for the very light axion can arise as a com- 
posite meson with an extra confining force whose scale is 
much larger than the electroweak scale. This confining 
force can be the hidden sector gauge group in supcrgrav- 
ity or just an extra gauge group. Let us call this extra 
confining gauge group 'axi-color' SU(AT). To create the 
QCD axion below the axi-color scale, there must be two 
classically conserved axial global symmetries |9l], 12261 
Kim (1985), Choi, Kim (1985)]. With only one axial 
symmetry, there would not result a masslcss meson, as 
in the case of one flavor QCD there is no massless meson 
since the only meson rj' becomes heavy b y th e instan- 
ton solution of the so-called U(l) problem J348I . 't Hooft 
(1976)]. For two axial symmetries, we can consider two 
kinds of axi-quarks, Q Aa , Qaoh 1 A i anc ^ Qa where A is the 
SU(iV) index and a is the SU(3) C index. For these vector- 
like representations, (N, 3) + (N, 3) + (N, 1) + (N, 1) un- 
der SU(7V)xSU(3) c , mass terms are not introduced. The 
axi-color vacuum angle problem is solved basically by the 
masslcss axi-quarks, Q and q. Even though Q looks like 
a masslcss QCD quark, it cannot be considered as the 
masslcss quark solution of the QCD 6 problem. After 
integrating out the axi-color degrees, we obtain an effec- 
tive Lagrangian resulting from Q and q. The axi-baryons 
are expected to be removed at the axi-color scale. Out of 
two kinds of mesons, one (the axi-color 77') is removed at 
the axi-color scale and the other remains exactly mass- 
less. However, this massless axi-color meson couples to 
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the QCD anomaly and becomes a QCD axion through 
the C3 term, becoming the so-called hadronic axion. Out 
of the two currents, 






(113) 



the divergence of J^ corresponds to the massless meson 
a below the axi-color scale, 
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and hence we obtain the effective interaction (11121). In 
this minimal model, the domain wall number is N J9ll.l92l. 
Choi, Kim (1985)]. In a supergravity model of preons, a 
similar mechanism was used to realize a composite axion 
[U, Babu, Choi, Pati, Zhang (1994)], where the role of 
q type matter is replaced by the meta-color gluino A' 
where meta-color is the binding force of preons. Even if 
the meta-color gluino obtains a mass of 0(100 GcV), the 
QCD 9 can be made within the experimental bound if F a 
is greater than 10 11 GeV. 

The composite axion of [lol, Chun, Kim, Nilles (1992)] 
is a composite made of hidden-color scalars whose bilin- 
ears develop VEVs and break the PQ sym metry. This 
idea has been made more concrete in [242|, Kim, Nilles 
(2009)]. 

In the gauge mediated SUSY brea king scenario a la 
Intrilligator, Sciberg and Shih (ISS) J204I . Intrilligator, 
Seiberg, Shih (2006)], for example, an SU(iV c ) confining 
group with Nf flavors satisfying N c + 1 < Nf < 
allows a SUSY breaking local minimum. If Nf 
3 (for example N c = 7 and Nf — 10) with one type 
Q Aa + Qau an( l Nf — 3 flavors of the type q A + q A , 
then there can exists a suitable local minimum where the 
composite axion envisioned in Eq. Q113P can be realized. 
In this case, the SUSY breaking scale and the composite 
axion scale are related as first tried in J225I . Kim (1984)]. 
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Queve do (1 999)]. With the Randall-Sundrum type warp 
factor [3131 . Randall, Sundrum (1999)], it is possible to 
introduce the intermediate scale with a Planck scale M F 
via t he Gi ddings-Kachru-Polchinski stabilization mecha- 
nism [172, Giddings, Kachru, Polchinski (2002)]. 

Here, we look only at the possibility of TeV scale M F . 
Since the Planck mass is given by M P « M F (RM F ) n / 2 , 
we obtain n > 2 for M F ~ 10 TeV QM HI Hannestad, 
Raffelt (2002), Kanti (2008)]. The Lagrangian in (4 + 
n)— dimensions ((4 + n)D) with a bulk field axion can be 
written as [75|, Chang, Tczawa, Yamaguchi (2000)], 

C cS = J d n yi \M%[d^adf*a + d y ad v a] 
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where vpq is the PQ symmetry breaking scale at the 
fundamental scale order M F , and 
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The 4-dimensional PQ symmetry breaking scale is F a « 

(Mp/M F ) A / n v PQ where A = |n| = ^Jn\ + Vn 2 n < n 

and F a falls between vpq and Mp. The very light axion 
is the n = component in Eq. (|116p . and the rest are the 
KK axions. The mass splitting of the KK axions are of 
order 1/-R, and the phenomenology of th ese K K axions 
for clkk ~^ 27 has been nicely studied in |l3li Di Leila, 
Pilaftsis, Raffelt, Zioutas (2000)], from which we have 
1/R - 1(10) cV for n = 2(3) for M F w 1 TeV. 

The possibility of a Z2 odd 5D gauge field in a warped 
fifth dimension has been suggested for a QCD axion 
under the assumption that all unwanted PQ symmetry 
breaking effects are suppressed [83|, Choi (2004)]. One 
such constraint is that the bulk fields carry the vanishing 
PQ charge. 



C. Axions with extra dimensions 

With large extra dimensions, the axion identification 
involves a few parameters: the fundamental scale mass 
M F , the Kaluza-Klcin (KK) radius R and the number of 
extra dimensions n. In addition, there are several ways 
to allocate the field(s) containing the axion in the bulk 
and/or branes. 

The possibility of large extra dimensions have been 
considered for the flat and warped extra dimensions. The 
TeV scale for M F was the main motivation to look for the 
next level of the current experimental limit on the mm 
scale gravity [12 . KL3l . Arkani-Hamed, Dimopoulos, Dvali 
(1998), Antoniadis, Arkani-Hamed, Dimopoulos, Dvali 
(1998)]. Because the axion scale is considered to be at 
the intermediate scale, a string theory at the intermediate 
scale M F has been considered also [54l . Burgess, Ibancz, 



D. SUSY breaking scale, axion and axino 

The 4D supergravity interactions with t he va nishing 
cosmological constant was written in 1983 |l!4l Crem- 
mer, Ferrara, Girardello, van Proyen (1983)]. The PQ 
symm etry can be embedded in the supergravity frame- 
work [Hi, Kim (1984)], 



2\ryl 



Wpq = (f 5 A 1 A 2 - Ff)Z_ + {f e A 1 A 2 - F*)Z- 
+.f Q A 1 Q 1 Q 2 
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where Z, Z', A\ and A2 are gauge singlet chiral fields, 
Q 1 and Q2 are chiral quark superfields, fs,f e ,Fi, and 
F$ are parameters. The superpotential (|117[) leads to 
the F-term SUSY breaking and the PQ symmetry break- 
ing at a common scale at order C(Ff,F|) if fg/f e 7^ 
Fi/F-2- The fq term defines the PQ charge of the 
heavy quark and the resulting axion is the KSVZ type. 
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The PQ symmetry breaking scale is given by nonzero 
(A 1 A 2 ) ~ (F 2 /\)cos([3 - a) and the SUSY breaking 
scale is given by nonzero Zp = —F 2 sin a sin (/3 



o 



^ 



F sinasin(/3 — a) where A 



and 

VJJTJ!,f 2 = 

X /-F 1 4 + F 2 4 ,tana = /e//« and tan/3 = Ff/Ff The ax- 
ino docs not obtain mass at this level, but obta i ns a 
mass at order the soft SUSY breaking scale [l02l . ll04L 
Chun, Kim, Miles (1992), Chun, Lukas (1995)]. Note 
that (Z'/Z) ~ -cot a with (Z,Z') = Q(F 2 JM). The 
early discussions on axino can be found in [1581 . Frcrc, 
Gerard (1983)]. 



E. The fi problem 

If Higgs doublets H U: d carry vanishing U(l) charges be- 
yond the MSSM gauge charges, then the superpotential 
can contain W^ = —fiH u Hd term where \x can be of order 
the fundamental scale since it is a supersymmetric term. 
This is problematic for the TeV scale electroweak sym- 
metry breaking, which is the so-called fj, problem [239l 
Kim, Nillcs (1984)]. This \i term is a supersymmetric 
Higgsino mass term and can be forbidden in W by in- 
troducing some symmetry, continuous or discrete. The 
widely discussed ones are the PQ and R symmetries. In 
the supergravity framework, if the Higgs doublets carry 
one unit of PQ charge then the nonrenormalizable inter- 
actions of the form S 2 H u Hd/Mp can be present in W if 
S 2 and H u Hd carry the opposite PQ charges. Then, the 
resulting p, is of order F 2 /Mp which can be of order of the 
TeV scale [HI, Kim, Miles (1984)]. With an intermedi- 
ate hidden sector with hidden sector squarks Q\ and Q 2 , 
one may have a nonrenormalizable interaction of the form 
QiQ 2 H u Hd/Mp. In this case also, the hidden sector 
squark condensation at the intermediate mass scale can 
generate a TeV scale \i [Tol, Chun, Kim, Miles (1992)]. 
[For the B^ term, one may consider {Q\Q 2 lMp}H u Hd 
in the Kahler potential.] The superpotential is better 
to poss ess this k i nd of PQ symmetry and/or R symme- 
try [63, [HE EM [H3, Hf Dine, Maclntire (1992), Hall, 
Randall (1991), Casas, Munoz (1993), Kim, Miles (1994), 
Kim(1999a)]. If so, even if the nonrenormalizable inter- 
actions are not considered, the gravity mediation scenario 
can ge nerat e a TeV scale /i via the Giudicc-Masiero mech- 
anism [l77l Giudice, Masiero (1988)]. In supergravity, 
the Higgsino mass ter m is prese nt in the chiral fermion 
mass matrix given by [114 , 12831 . Cremmer, Ferrara, Gi- 
rardello, van Proven (1983), Miles (1984)] 

e' G [G ,J - && - G^G-^fG^xuXLj (118) 

where G = K{^<j>*) - \n\W\ 2 . The term e-°G^ gives 
H ~ m 3 / 2 if K contains H u Hd [l77l Giudice, Masiero 
(1988)] and (j, ~ S 2 /M P if W contains S 2 H u H d /M P [HJ 
Kim, Miles (1984)]. 

In the ncxt-to-MSSM(NMMSM) models with W = 
SH u H d , the fi term can b e genera ted by the singlet 
VEV at the electroweak scale [67ll264l . Cerdeho, Hugonic, 



Lopez- Fogliani, Muhoz, Teixeira (2004), Lopez- Fogliani, 
Munoz (2006)]. In a Z'-added MSSM(Z'MSSM), the fi 
term can be successfully generated also [254L Langacker, 
Paz, Wang, Yavin (2008)]. 

Extending the MSSM gauge group which can be bro- 
ken down to the MSSM at a high energy scale, one can 
generate a reasonable \x. For example, there exists an in- 
teresting solution to the problem, "Why is there only one 
pair of Higgsino dou blets at low energy?" , in the extended 
electroweak model [257J Lee, Weinberg (1977)]. This is 
dictated by the extended gauge symmetry. This one pair 
problem is elegantly solved in the SUSY Lee- Weinberg 
type mod el due to the antisymmetric Higgsino mass ma- 
trix J234I . Kim (2007)], reminiscent of the "color" intro- 
duction to put the low-lying baryons in the completely 
symm etric representation 56 in the old flavor-spin SU(6) 
(l87l Han, Nambu (1965)]. 

Thus, explicit models toward a successful fi in the grav- 
ity mediation scenario can be constructed in extra singlet 
models, in SUSY-GUT models, through the superpoten- 
tial, through the Kahler potential, and in composite mod- 
els. 

The loop effects are the important sources of SUSY 
breaking in the gauge mediated SU SY breaking scenario 
(GMSB) [TH, Eslllll EH EM Dine, Nelson (1993), 
Dine, Nelson, Shirman (1995), Dine, Fischler, Srednicki 
(1981), Dimopoulos, Raby (1981), Dine, Fischler_ 
in anomaly mediation SUSY breaking (AMSB) 
Randall, Sundrum (1999), Giudice, Luty, Murayama, 
Rattazzi ( 1999 )], and even in the mirage mediation sce- 
nario [84|, 12651 Choi, Jcong, Okumura (2005), Loaiza- 
Brito, Martin, Miles, Ratz (2005)]. The GMSB has 
been suggested to solve the flavor problem [161] , Gab- 
biani, Gabrielli, Masiero, Silvestrini (1996)], present in 
the gravity mediation scenario. In this GMSB or any 
other loop generated SUSY breaking scenarios, the soft 
terms generated by supergravity effect are required to be 
sub-dominant compared to those arising from the loops 
or at best comparable to it. If it is sub-dominant as in 
the GMSB or AMSB, then there are some problems. 

Firstly, the generation of \x is difficult because p\, term 
generation via the Giudice-Masiero mechanism is sub- 
dominant at the TeV scale. One has to generate /x em- 
ploying the PQ and/or R symmetries, which however 
does not belong to the grand design of generating all TeV 
scale parameters dynamically. In this regard, another 
confining group around TeV scale has been proposed [88l . 
Choi, Kim (2000)], and the model presented there be- 
longs to a composite SU(2)w axion discussed in Subsec. 
IVIB1 which w as sa ved by introducing singlets and rele- 
vant couplings 266, Luty, Tcrning, Grant (2001)]. Then, 
again it does not succeed in generating all TeV scale pa- 
rameters dynamically. 

Second, in the loop SUSY breaking scenarios for gen- 
erating all TeV scale electrowe ak p arameters by loops, 
there exists the B^j ' \i problem [1521 . Dvali, Giudice, Po- 
marol (1996)]. Since it occurs at loop orders, we consider 
fd 4 6H u H d X^ for p, and fd 4 8H u H d XX^ for B^ where 
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the auxiliary component of X develops a VEV. From 
this observation, one generically obtains S„ ~ /iA where 
A ~ A*// 2 can be greater than /i which w as re medied 
by making B^ appear at a two loop order [1521 . Dvali, 
Giudicc, Pomarol (1996)]. This B^/fi problem occurs 
basically from the difference of the engineering dimen- 
sions of the Bfj, and ii terms. Both generically appear at 
one loop order with the coefficient g 2 /16ir 2 , and hence 
in describing the electroweak scale B^ term lacks one 
power of g 2 /16n 2 . Recently, a better solution employing 
a Kahl er p otential H u H,i{k\X + InX^) has been sug- 
gested |l75l . Giudice, Kim, Rattazzi (2008)], which can 
be compared to the original Giudice-Masiero Kahler po- 
tential H u HdX^ + ■ ■ ■ . There exist several m ore ideas to- 
ward the B^/n problem [zl, [To3, Hzl HI Cohen, Roy, 
Schmaltz (2007), Murayama, Nomura, Poland (2008), 
Roy, Schmaltz (2008), Cho (2008)] 

Maybe, the nonrenormalizable interactions are the 
easy solution of the (i and B^j \i problems even in the 
GMSB also. Here however one introduces another scale. 
Without a detail knowledge on the ultra-violet comple- 
tion of the MSSM, the nonrenormalizable interactions are 
usually assumed to be suppressed by the Planck mass 
Mp. But, there can be some heavy mass scale M, which 
can be somewhat smaller than the Planck mass Mp, for 
the seesaw mass of the required nonrenormalizable inter- 
actions. In string compactifications, it is known explic- 
itly that M can be different from M P [M S HH Choi, 
Kim (2006), Kim, Kim, Kyae (2007), Choi, Kobayashi 
(2008)]. A simple diagram giving an M dependence is 
shown in Fig. |2"61 where SU(2) doublets Hi and H 2 forms 
a vectorlike superheavy pair. This is a kind of the seesaw 



e.g. through the PQ symmetry preserving term in the 
Kahler potential 
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FIG. 26: The generation of the fj, term by a seesaw mecha- 
nism. 

mechanism of Higgsino doublet pairs. For this scenario, a 
superpotential possessing the PQ symmetry can be con- 
structed 
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where (X) = M and H u Hd term is forbidden by the PQ 
symmetry. Then, the ii term for H u and Hd (which give 
mass to up and down type quarks, respectively) is given 
by H = (fif 2 S 2 /M) [23l, Kim, Nilles (1984)]. If (S) 
is lowered to the hidden sector confining scale of order 
~ 10 10 ~ 12 GeV in the GMSB, the Higgsino mass can be 
made around the TeV scale by adjusting fif 2 /M. One 
may construct models with appropriate F-terms such 
that J5 M and m 2 oft are of the same order in the GMSB, 
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which also gives a li term. 

In the so-called "mixed mediation(M-mediation) sce- 
nario" , with comparable moduli, anomaly and gauge me- 
diations, that includes in its parameter space the GMSB, 
the AMSB, the mirag e mediation, and the deflected 
mirage mediation [1571 . Everett, Kim, Ouyang, Zurck 
(2008)], the loop generated \x term in general has a se- 
vere J5„/ ' ii problem. It seems that the model presented in 
an AMSB scenario [HI HI Pomarol, Rattazzi (1999), 
Rattazzi, Strumia, Wells (2000)] has the basic ingredient 
for the solution of B g (jx problem according to a PQ sym- 
metry as stressed in |l75l Giudice, Kim, Rattazzi (2008)]. 
This can be gleaned also from the axio n shift symmetry 
in the mirage mediation scenario [280J, Nakamura, Oku- 
mura, Yamaguchi (2008)]. 



F. Axions from superstring 

The most interesting theory housing axions is super- 
string theory Axions from string are described by ef- 
fective field theory below the compactification scale. If 
it arises from the spontaneous symmetry breaking of a 
tree-level global symmetry as we discussed so far, the 
answer is simple: There is no such axion since string the- 
ory would not allow any global symmetry If the com- 
pactification process leads to the SM, the renormalizablc 
terms in this effective theory respect the gauge symmetry 
SU(3) c xSU(2)xU(l)y and the global symmetries: the 
baryion number U(l)s and the separate lepton numbers 
U(1)lj. On the other hand, if the nonrenormalizable 
terms are allowed, one can write for example QlIlurcIr, 
breaking both the baryon number and the lepton number. 
Including the nonrenormalizable terms, thus the SM does 
not respect the baryon and lepton numbers. Similarly, 
there is no PQ global symmetry if we are allowed to write 
all nonrenormalizable terms. For the PQ symmetry, the 
situation is more severe. Suppose the singlet carrying 
the PQ charge is a. Then, a* a respects the PQ sym- 
metry but a 2 and a* 2 do not respect the PQ symmetry, 
which has led to the discussion on the gravitational ef- 
fects on the axion [H EH, GlS EM HD, Kamionkowski, 
March-Russell (1992), Barr, Seckel (1992), Holman, Hsu, 
Kephart, Kolb, Watkins, Widrow (1992), Ghigna, Lusig- 
noli, Roncadelli (1992), Dobrescu (1997)]. Therefore, the 
PQ symmetry cannot be discussed in general terms in 
terms of matter fields only, when we include gravity in 
the discussion such as in string theory. 

Thus, in string compactification one must consider the 
gravity multiplet also. Here, the gauge singlet bosonic de- 
grees in the gravity multiplet are graviton gMN (M, N = 
0, !,■•■ ,9), antisymmetric tensor Bmn and dilaton <J>. 
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In 10D, gMN and Bmn are gauge fields. A 4D action on 
Minkowski space x^(fi = 0,1,2,3) is obtained by com- 
pactifying six internal spaces y % (i — 4, • • ■ ,9) with the 
compact volume Vz- Some bosonic degrees from the 10D 
antisymmetric tensor field behave like pseudoscalars in 
the 4D effective theory. Thus, the axion candidates, if 
they are not arising from the matter multiplets, must be 
in Bmn- The pscudoscalar fields in Bmn arc like phase 
fields in axion models in field theory. Because there is 
no global symmetry in string theory, there must be no 
massless Bmn, otherwise the shift symmetry of Bmn 
must have worked as a global symmetry. From the tree- 
level equations of motion all pseudoscalar Bmn fields are 
not massless. For example, if a shift symmetry of Bmn 
is related to an anomaly as in the PQ current case, we 
consider that the shift symmetry is already broken. In 
other words, there is no shift symmetry of pscudoscalar 
Bmn if it is not anomalous. 

One must deal with these bosonic degrees in string 
compactification to see whether these components do 
not lead to terms in the potential (or the superpotential 
in SUSY models), which is very technical and model- 
dependent procedure. Here, we briefly discuss the story 
on axions from string and comment on its phenomeno- 
logical viability. Some r elevant r ecent reviews describing 
details can be found in flQ8ll343l. Svrcek, Witten (2006), 
Conlo n (20 06)1. The M-theory discussion was presented 
in [HIM! Choi (1997), Svrcek, Witten (2006)]. 

The pscudoscalar fields in Bmn come in two cate- 
gories, one the tangential component B^ v and the other 
Bij. B^v can be discussed in any string compactifica- 
tion and hence is called 'mo del- independent '(MI) while 



B^ depends on the compactification scheme as its in- 
ternal coordinates i and j imply and is called 'model- 
dependent'(MD). After presenting the string formulae 
containing Bmn, we will discuss the MI axion in Sub- 
subsec. IVIF II and then the MD axions present in much 
more speculative models in IVIF 21 One can go directly 
to Subsubsec. IVIF II if not interested in the details on 
the origin of the couplings. 

Now , the re exists the standard formulae for the string 
action [29Sl Polchinski (1988), Eq. (13.3.22JL which was 
lacking in the early days of string axions [83, 13721 . Wit- 
ten (1984), Choi, Kim (1985)]. The Type II dilaton (f> u 
and coupling gu = e^" are related to the 10D gravita- 
tional coupling Kio by Mf = 1/k 2 = \k j g 2 u t\ where 
a' = £ 2 J(2n) 2 . For the Type II string, there are NS- 
NS and R-R fluxes which can give anomalous couplings. 
These co mplic ated system housing pseudoscalars is re- 
viewed in jl08l . Conlon (2006)] with a possible conclusion 
of the difficulty of realizing a QCD axion i n str ing mod- 
els with a workable moduli stabilization [210i . Kachru, 
Kallosh, Linde, Trivedi (2003)]. In heterotic string mod- 
els, there does not exist a reasonable moduli stabilization 
mechanism, even though an ambitious attempt has been 
proposed [33|, Becker, Becker, Fu, Tseng, Yau (2006)]. 
We however discuss heterotic string axions below because 
string axions were found first in heterotic string mod- 
els and key couplings on axion phenomenology might be 
similarly discussed also for the Type II string. As in the 
Type II string, the heterotic string coupling is related 

to the dilaton $ as gt = e®. The kinetic energy terms 

of gMN, Bmn and Am are [2981 . Polchinski (1988), Eqs. 
(12.1.39) and (12.3.36)], 
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where Tr v is the trace over vector representation and the 
Chcrn-Simons 3-form is 



U>3 
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For EsxEs, there is the adjoint representation and we 
use ^jTr a in place of Tr v . For the compact internal vol- 
ume Vz, the Planck mass is Mp = inVz/g 2 ^ and the 4D 
gauge coupling constant is g YM = ^ 7r g 2 ^/Vz or olym = 
g 2 £®/Vz- In most compactifications, the SM gauge fields 
arise from the level k = \ embedding and the coupling 
olym is the coupling strength at the compactification 
scale. If the SM gauge fields are embedded in the level 



k, the SM gauge coupling at the compactification scale 
will be smaller by the factor k. For interactions of Bmn, 
we consider the Bianchi identity, the gauge invariant cou- 
plings of gaugino \ |l27lll4Ct Derendinger, Ibanez, Nillcs 
(1985), Dine, Rohm, Seiberg, Witten (1985)], and the 
Green-Schwarz terms |l8l[ Green, Schwarz (1984)] 



dH 



16 



y(trRf\R-trFf\F), H M npxT mnp X , 

BAtrf AFAfAF+--- (123) 



where Hmnp is the field strength of Bmn , F is the field 
strength of the gauge field A, the gauge invariant fcrmion 
coupling is the SUSY counterpart of the relevant terms of 
Eq. (|12ip . and • • • denotes more Green-Schwarz terms. 
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It was argued that the Hmn p coupling to the gaugino 
must be a perfect square [l40l . Dine, Rohm, Seiberg, Wit- 
ten (1985)], which gives a vanishing cosmological con- 
stant even for a nonvanishing gaugino condensation with 
nonzero (Hmnp) |l27l . Il40l . Derendinger, Ibanez, Nillcs 
(1985), Dine, Rohm, Seiberg, Witten (1985)]. 



1. Model-independent axion 

B^y with [i and v tangent to the 4D Minkowski space- 
time is th e MI axion present in all string compactifica- 
tions [37lL Witten (1984)]. Because it is a 4D gauge bo- 
son, one cannot write potential terms in terms of B^ y 
and it is massless if one neglects the anomaly term. 
The number of transverse degrees in B^ v is one and 
can be expressed as a pseudoscalar a by dualizing it, 
H^yp ex F a €p l y pa d' J a. Even though it is massless at this 
level, the Bianchi identity of (|123[) gives an equation of 
motion of a as d 2 a = (l/32Tr 2 F a )G a ^G atny which hints 
that a might be an axion. For it to be really a QCD 
axion, ci + C3 should be nonze ro as discussed in Subsec. 
ImBl It is known that c 3 = 1 Hzl, Witten (1985)] with 
C3 defined in Eq. (flT)|) . The other possible couplings are 
given by the second term of (|123[) 
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which is the c\ term defined in Eq. (|19p . There is no 
ci term, and ci + c-s — 1 and hence -ff M „ p is really an 
axion and is model-independent. This is a hadronic ax- 
ion. This MI hadronic axion can have a nonvanishing 
c\ and hence its phenomenology can be different from 
that of the KSVZ hadronic axion. So, in Eqs. (|61l62p 
for the MI hadronic axion, one has to add the relevant 
Ci term from Eq. (|124p . The domain wall number of 
the MI axion has been shown to be Ndw = 1 by con- 
sidering the coupling of the MI axion to stri ng X M (a, r) 
on the world sheet / d 2 ae at3 B^d a X^dpX v [Si Witten 
(1985)]. F a is about 10 -3 times the Planck mass [89|, 
Choi, Kim (1985a)], and the correct relation, obtained 
from (fT2T|) . is F a /k = a c M P /2 3 / 2 n ~ 10 16 GeV where 
k is the level of t he SM embedding and a c is the QCD 
coupling constant |343l . Svrcek, Witten (2006)]. But the 
value F a ~ 10 16 GeV most probably overdoses the uni- 
verse. 

An idea for lowering the MI axion decay constant may 
be the following. In some compactification schemes, an 
anomalous U(l) an gauge symmetry results, where the 
U(l)an gauge boson eats the MI axion so that the U(l) a n 
gauge boson becomes heavy. This applies to the MI 
axion since the coupling d^dMiA*" 1 is present by the 



Grecn-Schwarz term [lol HlO, Witten (1984), Chun, 
Kim, Nilles (1992)]. In fact, even before considering 
this anomalous U(l) a ngauge boson, the possibility was 
pointed out by Barr [29|, Barr (1985)], which became a 
consistent theory after discoveri ng the ano malous U(l) an 
from string compactification [la . I136L [l43L Dine, Seiberg, 
Witten (1987), Atick, Sen (1987), Dine, Ichinose, Seiberg 
(1987)]. Then, a global symmetry survives down the 
anomalous U(l) an gauge boson scale. A detailed scenario 
is the following. The way anomalous U(l) an with the 
gauge transformation, # an — > # an +constant, is obtained is 
by calculating U(l) an charges of fermions. Thus, we have 
a non- vanishing 02 in Eq. (fl"9)) as rm/j L ipR exp(ic2# an ) and 
C3 of the MI axion as c^Omi {FF}. For all gauge group 
factors, the anomaly unit s ar e calculated and they are 
shown to be identical [63|, |228| . Casas, Kathehou, Munoz 
(1988), Kim (1988)]. For the MI axion to be the part of 
a gauge boson, it must be a true Goldstone boson with- 
out an anomaly, i.e. it should be exactly massless; so let 
us transform away the C3 term by a phase redefinition of 
fermions such that c-i = c<i — Cz(8mi) /$an and C3 = can 
be made for all gauge fields, i.e. ami coupling to anoma- 
lies vanishes for all gauge groups. Because the longitudi- 
nal gauge boson aui is removed away, we are left with 
the C2 term only, mip L ipuexp(ic20 an ) + h.c, without the 
need to consider the gauge symmetry U(l) an . At low en- 
ergy, however, the term m?/>£'0-R e xp(*C2# a n) has a global 
symmetry, f? an — > # an + constant, with # an not depending 
on x^ . Thus, the interaction mip L tpRexp(ic20a.n) + h.c. 
explicitly shows a global U(l) axial symmetry or the 
PQ symmetry below the U(l) an gauge boson mass scale: 



l-l'L 



^e-^n/2 and q p R _, ^e^U/2. T h is global P q 



4 Never thele ss, its property may depend on models in warped 
space Ill9l . Dasgupta, Firouzjahi, Gwyn (2008)]. 



symmetry can be broken in the axion window as in the 
field theoretic axion models. However, this idea on the 
decay constant does not work necessarily because most 
fields, including those removed at the GUT scale, carry 
the U(l)an charge. 



2. Model- dependent axion 



In 4D, Bmn contains more pseudoscalars Bij with i 
and j tangent to the compact space Vz ■ If they are ax- 
ions, these are MD axions. The number of massless B^ 
modes at the KK mass lev el is the second Betti num- 
ber of the compact space |l82l . Green, Schwarz, Wit- 
ten (1987), Eq. (14.3. 10)1, which was discussed in the 
early days in [90|, HzE Wa Witten (1984,1985), Choi, 
Kim (1985)]. The string propagation on M4 x Vz can be 
described by a suitable nonlinear sigma model. In this 
sigma model description, when a closed string wraps Vz 
topologically nontrivially then there are world-sheet in- 
stantons due to the map S\ — ► £7(1). It was known that 
the world-sheet instantons are pres ent precisely if the sec- 
ond Betti number is nonzero [1821 . Green, Schwarz, Wit- 
ten (1987)], and hence the MD axions are expec ted to re- 
ceive non-negligible masses non-perturbatively |l4ll . Il42l 
iH, Wen, Witten (1986), Dine, Seiberg, Wen, Witten 
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(198 6,1987)1, but it may be a model dependent statement 
[299l Polchinski (2006)]. If a MD axion is known to have 
no potential term except the anomaly terms, then one 
should check the c 2 and C3 couplings to confirm that it is 
really an axion. There has been no example presented yet 
in this way for a MD axion. If a MD axion is present, its 
decay constant is expected to be near the string scale as 
explicitly given as Fmd = ctj Mp/2 3 / 2 Trk 1 / 3 g s ' from 
the anomaly term alone in [3431 . Svrcek, Witten (2006)]. 
The Green-Schwarz term integrated over Vz leads to this 
kind of the decay constant for the MD axion [90, Choi, 
Kim (1985)]. However, as commented above, one has to 
calculate the corresponding c 2 term also to pinpoint the 
MD axion decay constant Fmd- 



3. Toward a plausible QCD axion from string 



A key problem in string axion models is to find a 
method obtaining a QCD axion at the axion window 
(10 9 GeV < F a < 10 12 GeV) but an attractive model 
toward this direction is still lacking. Thus, the most 
pressing issues is the problem of introducing a detectable 
QCD axion from superstring. It includes the search for 
an approximate PQ symmetry toward a detectable QCD 
axion. 

The conditions for compactificd manifolds in warped 
space toward low ering the MI axion decay constant has 
been discussed in [119i . Dasgupta, Firouzjahi, Gwyn], but 
its realization seems nontrivial. 

The idea of localizing MD axions at fixed points towa rd 
lowering the decay constant has been proposed in [2231 . I. 
W. Kim, J. E. Kim (2006)]. It uses the warp factor idea 
and on e ne eds a so-called Giddings-Kachru-Polchinski 
throat [ml, Giddings, Kachru, Polchinski (2002)] in the 
Type-II string, but in the heterotic string a non-Kahler 
Vz is needed [33, Becker, Becker, Fu, Tseng, Yau (2006)]. 
Indeed the warp factor is obtained in this way, which 
however has the power law behavior. 

The intermediate scale string models can introduce the 
axion window as the ultra-violet completion scale [54], 
Burgess, Ibahez, Quevedo (1999)]. On the other hand, in 
this case the large radius to generate the Planck mass is 
the scale needing explanation. 

We note that an attractive method of obtaining F a in 
the axion window is through the composite axion from 
superstring as discussed in Subsec. IVIBI However, from 
string construction the composite axion has not been ob- 
tained so far. 

Even if one lowered some F a , we must consider the 
hidden sector also in estimating the axion masses and 
decay constants as discussed below. 



4- Hidden sector confining forces, axion mixing, and 
approximate PQ symmetry 

With the hidden-sector confining forces, we need at 
least two (QCD and one hidden-sector) 8s which have to 
be settled to zero and hence we need at least two axions. 
For definitcness, let us consider only one more confining 
force at an intermediate scale, which may be the source 
of gravity mediation or GMSB. In this case, at least one 
MD axion is assumed to be present, and the axion mixing 
must be considered. We assume that one decay constant 
is in the intermediate scale. Here, there is an important 
(almost) theorem: the cross theorem on decay constants 
and condensation scales. Suppose that there are two ax- 
ions a\ with iq and a 2 with F 2 (iq <C F 2 ) which couple to 
axion potentials wit h scales Ai and A2 (Ai <C A2). The 

theorem states that [22IS3IHH Kim ( 1999 ' 2000 )' L 
W. Kim, J. E. Kim (2006)]: according to the diagonal- 
ization process in most cases with generic couplings, the 
larger potential scale A2 chooses the smaller decay con- 
stant iq, and the smaller potential scale Ai chooses the 
larger decay constant F 2 . So, just obtaining a small de- 
cay constant is not enough. The hidden sector may steal 
the smaller decay constant. Probably, the QCD axion 
is left with the larger decay constant. We can turn this 
around such that the hidden sector instanton potential is 
shallower than the QCD instanton potential since the in- 
stanton potential is proportional to the light quark mass 
as discussed in Subsec. IIIIBI If a hidden-sector quark 
mass is extremely small, then the QCD axion can win 
the smaller decay constant and the other axion is an ex- 
tremely light axion wh i ch ca n be used to fit the observed 
dark energy [HI HH HH Perlmutter et. ai. (1998), 
Riess et. ai. (1998), Komatsu et. al (WMAP Collabo- 
ration, 2008)] which is named as the quintessential axion 
[24l Kim, Nilles (2003)]. This can be easily realized if 
some hidden-sector squark condensations are very small 
as Fig. [57] can generate hidden-sector quark masses [2231 . 
I.-W. Kim, J. E. Kim (2006)]. 
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FIG. 27: The hidden-sector squark condensation breaks chiral 
symmetry and generate hidden-sector quark masses. 

Since it is difficult to obtain a reasonable light MD ax- 
ion, it has been tried to look for an approximate PQ sym- 
metry from string compactification. Only one reference 
exists from a realistic string compactification because of 
the difficulty of calculating all approximate PQ charges of 
quarks 0, Choi, I.-W. Kim, J. E. Kim (2007)]. After all, 
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the topologically attractive Bij may not be the QCD ax- 
ion we want. In this regard, we note that there already 
exists a field theoretic work regarding an approximate 
PQ symmetry, starting with a discrete Z$ symmetry 
J255I . Lazarides, Panagiotakopoulos, Shafi (1986)]. Later, 
the gravitational nonpcrturbative effects such as worm- 
holes and black holes were pheno menq l ogica l ly st udied 
in view of any global symmetries |l7li Il74l . l260l Gid- 
dings, Strominger (1988), Lee (1988), Gilbert (1989)]. It 
is known that the PQ symmetry breaking operators in 
the supe r potential must be forbidden up to dimension 8 
[H ELm 1l7Cl [Hi, [HI Kamionkowski, March-Russell 
(1992), Barr, Seckel (1992), Holman, Hsu, Kephart, 
Kolb, Watkins, Widrow (1992), Ghigna, Lusignoli, Ron- 
cadelli (1992), Dobrescu (1997)]. If we introduce an ap- 
proximate PQ symmetry, it is better to forbid the PQ 
symmetry breaking operators up to dimension 8 in the 
superpotential, possibly up to dimension 7 with reason- 
ably small couplings somewhere. 

In this spirit, it is worthwhile to check approximate 
PQ symmetri es in string derived models. The MSSMs 
presented in HH, HM ^M, HH Kim, Kyae (2007), J. 
E. Kim, J.-H. Kim, Kyae (2007), Kim (2007a, 2007b)] 
satisfy most phenomenological constraints and one can 
check approximate global symmetries. But it is a tedious 
work, and so far only for the model of |237l . Kim, Kyae 
(2007)] an approximate PQ symmetry has been checked 
out. In searching for an approximate global symmetry in 
a string derived model, there are so many Yukawa cou- 
plings to be considered that a complete study up to all 
orders is almost impossible. For example, in [94 , Choi, 
I.-W. Kim, J. E. Kim (2007)] there appear O(10 4 ) d = 7 
superpotential terms and it is not a trivial task to find 
an approximate PQ symmetry direction, considering all 
these terms. Up to dimension 7 terms, there exists an 
approximate PQ symmetry which is spontaneously bro- 
ken. The resulting axion coupling with photon has been 
calculated in 0, Choi, I.-W. Kim, J. E. Kim (2007)] 
which is presented in Fig. [16] to gether w ith the CAST 
and Tokyo axion search bounds [lOt 12031 . Andriamonje 
et. al. (2007), Inouc et. al. (2008)]. But the axion decay 
constant is not lowered. It is because the needed sin- 
glet VEVs, leading to the low energy MSSM, carry PQ 
charges. This is a generic problem for observable axions 
from superstring. Comparing this to the MI axion case 
with the anomalous U(l) an , it may be easier to realize the 
observable axion with an approximate PQ symmetry. 



VII. AXINO COSMOLOGY 

Supersymmetrization of axion models includes the 
fcrmionic superpartner axino a and the scalar superpart- 
ner saxion as discussed in Subsec. IVIDl Both saxion and 
axino masses are split from the almost vanishing axion 
mass if SUSY is broken. The precise value of the ax- 
ino mass depends on the model, specified by the SUSY 
breaking sector and the mediation sector to the axion 



supcrmultiplet [283, Nilles (1984)]. Most probably, the 
saxion mass is around the soft mass scale Mstjsy- The 
axino mass should also be near this scale as well. But 
the a xino mass can also be much smaller than .Mstjsy 
fiol OH, [HI, Frere, Gerard (1983), Kim, Masiero, 
Nanopoulos (1984), Chun, Kim, Nilles (1992a)] or much 
larger than M SUSY |l04 Chun, Lukas (1995)]. Therefore, 
we take the axino mass as a free parameter here. 

The decoupli ng t emperature of the axino supermulti- 
plet is of order [3 111 . Rajagopal, Turner, Wilczck (1991)], 
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10 11 GeV 



F a 



10 12 GeV 



£ ] d25) 
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where a c is the QCD coupling constant. 

The saxion cosmology is a simple extension of the 
standard cosmology with saxion mass around the SUSY 
breaking scale 0, |zl [HI J- E. Kim (1991), Chang, 
H. B. Kim (1996), Asaka, Yamaguchi (1999))], but 
its effect is not so dramatic as the effect of the ax- 
ino. Th e refo r e, he re we focus on the axino cosmology 
[H, OH, [HI, [311, Rajagopal, Turner, Wilczek (1991), 
Covi, J. E. Kim, Roszkowski (1999), Covi, H. B. Kim, J. 
E. Kim, Roszkowski (2001), Choi, J. E. Kim, Lee, Seto 
(2008)] . In the moduli stabilization scenario of Ref. [210l ] , 
the saxion VEV has been estimated in [85|, Choi, Jeong 
(2007)]. 

The axino cosmology depends crucially on the nature 
of R-parity. If R-parity is conserved and the axino is 
lighter than the neutralino, then most probably the axino 
or gravitino (in case of GMSB) is the LSP. If R-parity 
is not conserved, the neutralino can decay to ordinary 
SM particles, which has been discussed extensively in 
[8J, Allanach, Dedes, Dreiner (2004)] and references cited 
there. 

Now we focus on R-parity conservation. The neu- 
tralino, if it is the LSP, is a natural candidate for dark 
matter. Due to TeV scale sparticlc interactions, the ther- 
mal history of neutralinos allows them to be dark mat- 
ter. But, imposing a solution of the strong CP problem 
via the axion, the thermal history involves contributions 
from the axion sector, notably by the axino. Since axino 
cosmology depends on neutralino and gravitino number 
densities, let us comment on the neutralino and gravitino 
cosmologies before discussing the effect of the axino. The 
neutralino cosm ology dep ends on the neutralino freeze- 
out temterature [l48l 12581 . Lee, Weinberg (1977), Drees, 
Nojiri (1993)] and the gravitino cosmolo gy d epends on 
the reheating temperature after inflation [362, Weinberg 
(1982)]. Here we list several relevant temperatures in the 
axino cosmology 

Ta-dcp = axino decoupling temperature 

Xr = reheating temperature after inflation 

Tf r = neutralino freeze — out temperature (126) 



T„ 



rad 



axino — radiation equality temperature 



Td = radiation temperature right after a, decay 
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Here, we arc interested in the axino domination of the 
dark matter density. The cold axino dark matter possi- 
bility in the evolution history occurs either a heavy axino 
decayed already or axino has not decayed yet. If the ax- 
ino has not decayed yet, the current axino CDM can be 
estimated using Ta-dcp or Tr. If it has decayed already, 
the past cold axino dark matter requires the existence of 



;maYa(T™ n ) 



T r . 



(127) 



so that Ya{T R ) = n-a/s > Y a (T™ in ) at the time of the 
reheating after inflation, where T^ ln is the temperature 
above which axinos dominate the universe before they 
decay. 



to a stronger bound Tr < 10 8 GeV if the gravitino is 
lighter than the gluino and Tr < 10 7 GeV if the grav- 
itino is heavier than the gluino. This gravitino prob- 
lem is absent if the gravitino is the next LSP (NLSP), 



ma < m 3 / 2 < 



'v 



since a thermally produced gravitino 



would decay into an axino and an axion which would 
not affect the BBN produced light elements [la Asaka, 
Yanagida (2000)]. 

If the gravitino is the LSP with the stau or neutralino 
as the NLSP, the gravitino can be the CDM even in 
the constrained MSSM (or mSUGRA) for some param- 
eter space, avoiding the BBN and b — ► sj constraints 
IH HI, [l5l, [HI Boehm, Djouadi, Drees (2000), Ellis, 
Olive, Santoso, Spanos (2004), Roszkowski, Ruiz de Aus- 
tri, Choi (2005), Ccrdcno, Choi, Jedamzik, Roszkowski, 
Ruiz de Austri (2006)]. 



A. Neutralino and gravitino 



B. Axino 



The neutralino LSP seems the most attractive can- 
didate for CDM simply because the TeV or der SUSY 
breaking scale introduces the LSP as a WIMP [154 Il79l 
Goldberg (1983), Ellis, Hagelin Nanopoulos, Olive, Sred- 
nicki (1984)]. The neutralino, which was in the thermal 
equilibrium in the early universe, decouples and freezes 
out when the annihilation rate becomes smaller than 
the Hubble parameter. The fr eeze - out t emperature Tf r 



is normally given by m x /25 [246 



Lee, Weinberg 



(1977b), Kolb, Turner (1990)], e.g. 4 GeV for 100 GeV 
neutralino. Obviously, the neutralino relic density is not 
affected by the axino for Tb > Tf r since neutralinos were 
in thermal equilibrium after the axino decay. This is the 
standard neutralino dark matter. With the axino intro- 
duction, therefore, we study the case Trj < TV r . 

Gravitinos in the universe are important if they dom- 
inate the dark matter fraction now or they had affected 
the result of nucleosynthesis. Thermal gravitinos pro- 
duce d at the Planckian time are important if m 3 / 2 ~ 1 
keV HH Pagels, Primack (1982)]. However, in the in- 
flationary scenario these Planckian time gravitinos are 
not important now. For a heavy gravitino, it was ob- 
ser ved t hat the gravitino decay affects the nucleosynthe- 
sis |362l . Weinberg (198 2 )], w hich was suggested to be 
solved by inflation [HI, \2M, Krauss (1983), Khlopov, 
Linde (1984)]. Then, the gravitino number density is 
roughly estimated in terms of the reheating temperature 
after inflation n 3 / 2 oc Tr. To estimate the cosmological 
bound on T r rather accurately, a full supergravity inter- 
action [ll4l . Cremmer, Ferrara, Girardello, van Proyen 
(1983)] has been used and applied to the dissociation 
problem of rare light elemen ts su ch as deuterium, etc., 
resulting in T R < 10 9 GeV [HI, Ellis, Kim, Nanopou- 
los (1984)]. A recent calculation of Tr has been per- 
formed using the nucleosynthesis co de to look for 7 Li de- 
struction and/or 6 Li overproduction 217l . l220l . Kawasaki, 



Thus, in SUSY theories we must consider a relatively 
small reheating temperature 10 7-8 GeV. The axino cos- 
mology must also be considered with this low reheating 
temperature. 

In principle, the axion supcrmultiplet is independent 
from the observable sector in which case we may take 
the axino mass as a free parameter of from ke V scale to 
a value much larger than the gravitino mass [l02l Il04 
Chun, Kim, Miles (1992a), Chun, Lukas (1995)]. Light 
axinos (ma < 100 GeV) can be a dark matter can- 
didate, and they have been stud ied extensively as a 
warm dark matter candidate [3 111 Rajagopal, Turner, 
Wilczek (1991)] with the reheating temperature given in 



Kohri Moroi (2005, 2008)], following the earlier work 



of [116|, Cyburt, Ellis, Fields, Olive (2003)], which led 



51, Brandenburg, Steffen (2004)], or a CDM candidate 
H, [lOi, \lM, HM [Ml Covi, J. E. Kim, Roszkowski 
(1999), Covi, J. E. Kim, H. B. Kim, Roszkowski (2001), 
Roszkowski, Seto (2007), Seto, Yamaguchi (2007), Asaka, 
Yanagida (2000)]. Heavy axinos, however, cannot be the 
LSP and can decay to the LSP plus light particles. This 
heavy a xino decay to neutralinos has been already con- 
sidered |l04l Chun, Lukas (1995)]. The heavy axino pos- 
sibility was considered brief l y in studying cosmological 
effects of the saxion by |219l l22dl . Kawasaki, Nakayama, 
Senami (2007), Kawasaki, Nakayama (2008)]. A more 
complete cosmological analysis of the heavy axino has 
been discussed in [971 . Choi, Kim, Lee, Seto (2008)]. 

Si nce the CDM fraction of the universe is roughly 0.23 
[247l Komatsu et al. (WMAP Collaboration, 2008)], let 
us focus on the CDM possibility of the axino or axino 
related neutralino. 

For the axino to be the LSP, it must be lighter than 
the lightest neutralino and gravitino. In this case, we 
do not have Tb of Eq. Q126p . If the lightest neutralino 
is the NLSP, ma < m x < m 3 / 2 , the thermal produc- 
tion (TP) mechanism gives the afore-mentioned bound 
on the reheating temperature after inflation. At the high 
reheating te mpe rature, the TP is dominant in the axino 
production [lid Covi, J. E. Kim, Roszkowski (1999)]. 



42 




10 



10 -6 1Q -4 1Q -2 

m~ a [GeV] 



l() z 



FIG. 28: The solid line is the upper bound from TP. The 
yellow region is the region where NTP can give cosmologically 
interesting results (f2g TP /i 2 ~ 1). The freezeout temperature 
isT F «^. 
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FIG. 29: The T R vs. m a plot for m x = 100 GeV and F a 
10 11 GeV. 



If the reheating temperature is below the critical energy 
density line, there exists another axino CDM possibility 
from the nonthermally prod uced (NTP) axinos which re- 
sult from neutralino decay [l09l . Covi, J. E. Kim, H. B. 
Kim, Roszkowski (2001)]. This situation is sketched in 
Fig. [281 We note that with R-parity conservation, the 
double production of low mass axinos are negligible in su- 
pernovae, and hence from SN1987A there does not result 
a useful exclusion region in the low mass region. 

Since the final axino energy fraction is reduced by the 
mass ratio, Og/i = (rria/m x )fl x h' 2 for rria < m x < 
to 3/2j the stringent cosmologically constrained MSSM 
parameter space for m x can be expanded. As shown 
in this figure, the NTP axinos can be CDM for rel- 
atively low reheating temperature (< 10 TeV) for 10 
McV < iria < m x . In Fig. [28] the thin dash line cor- 
ner on the RHS corresponds to the MSSM models with 
Q, x h? < 10 4 , and a small axino mass renders the possibil- 
ity of the axino forming 23% of the closure density. If all 
SUSY mass parameters are below 1 TeV, then probably 
Q, x h 2 < 100 (the thick solid-dash line corner on the RHS) 
but a sufficient axino energy density requires mg > 1 
GeV. Thus, if the LHC docs not detect the neutralino 
needed for its closing of t he univer se, the axino clos- 
ing is a possibility [98|, [93, 1 1 111 . Ill2i Choi, Roszkowski, 
Ruiz de Austri (2008), Covi, Roszkowski, Small (2002), 
Covi, Roszkowski, Ruiz dc Austri, Small (2004), Choi, 
Roszkowski (2006)]. If the NLSP is stau with the axino 
or gravitino LSP, the previously forbidden stau LSP re- 
gion is erased. In this case, the CDM axino is similar to 
the bino LSP case but it is easier to detect the stau signal 
due to the charged stau at the LHC [5l|, Brandenburg, 



Covi, Hamaguchi, Roszkowski, Steffen (2004)] Ho wever, 
the efforts to detect axinos may be difficult [222l . H. B. 
Kim, J. E. Kim (2002)]. 

In the GMSB scenario, the gravitino mass is generally 
smaller than the neutralino mass and possibly smaller 
than the axino mass, for which cas e the cosmological ef- 
fect has been studied by [l0l|,[22l|, Chun, H. B. Kim, J. 
E. Kim (1993), H. B. Kim, J. E. Kim (1995)]. 

For a heavy axino decaying to a neutralino, we present 
a T R vs. m~ a plot for F a = 10 11 GeV in Fig. CSO The 
region T R > Ta_dcp is above the dashed blue line. The 
axino lifetime greater than 0.1 s is denoted by the red 
shaded region in the LHS. The blue shaded region in the 
RHS is where the axino decays before the neutralino de- 
couples (Td > Tf r ). The magenta lines (horizontal) are 
the contours of the entropy increase due to the axino 
decay, r = Sf/So. Above r = 1 lines axinos dominate 
the universe before they decay. The green lines (verti- 
cal) denote the (cr arm w re ;), where a ann is the neutralino 
annihilation cross section, in units of GeV - which are 
used to give the right amount of neutralino relic den- 
sity. In Fig. [29] we use neutralino and gluino masses 
as m x — 100 GeV and nig = 2 TeV, respectively. For a 
larger F a and a heavier neutralino mass, the green lines 
move to the right [97|, Choi, Kim, Lee, Seto (2008)]. 
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